arXiv:1509.05334vl [astro-ph.GA] 17 Sep 2015 


Stars, Gas, and Dark Matter in the Solar Neighborhood 
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ABSTRACT 

The surface density and vertical distribution of stars, stellar remnants, and gas in the 
solar vicinity form important ingredients for understanding the star formation history 
of the Galaxy as well as for inferring the local density of dark matter by using stellar 
kinematics to probe the gravitational potential. In this paper we review the literature for 
these baryonic components, reanalyze data, and provide tables of the surface densities 
and exponential scale heights of main sequence stars, giants, brown dwarfs, and stellar 
remnants. We also review three components of gas (H 2 , HI, and HII), give their surface 
densities at the solar circle, and discuss their vertical distribution. We find a local total 
surface density of M dwarfs of 17.3 ± 2.3 Mq pc“^, significantly higher than previous 
values. Our result for the total local surface density of visible stars (main sequence stars 
and giants), 27.0 ± 2.7 Mq pc“^, is close to previous estimates due to a cancellation 
of opposing effects: more mass in M dwarfs, less mass in the others. The total local 
surface density in white dwarfs is 4.9 ± 0.6 Mq pc“^; in brown dwarfs, it is ~ 1.2 
Mq pc“^, but with considerable uncertainty. We find that the total local surface 
density of stars and stellar remnants is 33.4 ± 3 Mq pc“^, somewhat less than previous 
estimates but within the errors of many of them. We analyze data on 21 cm emission 
and absorption and obtain good agreement with recent results on the local amount of 
neutral atomic hydrogen obtained with the Planck satellite. The local surface density 
of gas is 13.7 ± 1.6 Mq pc“^. The total baryonic mass surface density that we derive 
for the solar neighborhood is 47.1 ± 3.4 Mq pc“^ (43.8 Mq pc“^ within 1.1 kpc of 
the midplane). Combining these results with others’ measurements of the total surface 
density of matter within 1-1.1 kpc of the plane, we find that the local density of dark 
matter is pdm = 0.013 ± 0.003 Mq pc“^ = 0.49 ± 0.13 GeV cm“^. The local density 
of all matter is 0.097 ± 0.013 Mq pc“^. We discuss limitations on the properties of a 
possible thin disk of dark matter. 
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1. Introduction 

A central goal in the study of the origin and structure of galaxies is to understand how both 
baryonic and dark matter are distributed within them. For the Milky Way, essential input is 
provided by the volume density and surface density of stars, stellar remnants and gas in the solar 
vicinity—the local baryon budget. This input is required for accurately inferring the amount of 
dark matter in the solar vicinity by dynamical modeling of observations of stars in the gravitational 
potential of both the dark matter and the baryonic matter, whether this modeling is local, based on 
fitting observations of stars in the solar vicinity, or global, fitting the rotation curve of the Galaxy 
(Read 2014). The value of poM is important for estimating the sensitivity requirements of dark 
matter searches and for constraining the shape of the Milky Way dark matter halo, which can be 
affected by late major mergers (Read 2014). It has also been suggested that there is a thin disk 
(scale height h = 0(100 pc)~Shaviv et al. 2014) or very thin disk {h = 0(10 pc)-Fan et al. 2013; 
Randall & Reece 2014; Kramer & Randall 2015) of dark matter that brings the period of vertical 
oscillations of the Sun into agreement with the ~ 30 Myr period of different phenomena observed 
on Earth. 

The distribution of stars by mass (the present day mass function, or PDMF) provides important 
information on the star formation history of the Galaxy, and the distribution of stars by distance 
from the Galactic plane provides information on the dynamical evolution of the stellar disk. The 
shape of the PDMF also provides constraints on the stellar initial mass function (IMF). In this 
work, we shall assume that the IMF of stars is universal in galactic disks when averaged over large 
regions of space and long intervals of time (Bastian et al. 2010). Although significant variations 
in the stellar mass function of star clusters are observed (Dib 2014; see also Parravano, McKee, & 
Hollenbach 2011, hereafter PMHll), the summed mass distribution of objects from various clusters 
converge to a smooth shape, the universal IMF. The stars in the solar neighborhood have been 
formed in a huge number of events and environments during the lifetime of the local disk (e.g., 
Bland-Hawthorn et al, 2010), so cluster-to-cluster variations should not affect our analysis. No 
general trends of the shape of the mass function of stars in clusters on the size of the cluster or its 
metallicity have been clearly established. The only deviations from a universal IMF that have been 
established are in massive elliptical galaxies (e.g., Conroy et al. 2013, but see Smith et al. 2015) 
and near the Galactic Center (e.g., Lu et al. 2013), although it should be noted that there is no 
evidence for deviations from a universal IMF in the old nuclear star cluster at the Galactic Center 
(Fritz et al. 2014). 

A standard reference for the baryon budget in the solar neighborhood (which we define as a 
cylinder of radius 1 kpc extending above and below the plane) is Flynn et al. (2006; hereafter F06). 
We make three improvements on their results: First, we alter their PDMF so that it is consistent 
with the observations of stars within 25 pc of the Sun as determined by Reid et al. (2002; hereafter 
RGH02). This involves both increasing the surface density of M dwarfs based on a revision of the 
analysis of Zheng et al. (2001; hereafter ZOl) and decreasing the surface densities of more massive 
stars. Second, we determine the surface density of white dwarfs in two ways; We implement the 
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correction suggested by Katz et al (2014), who pointed out that a significant number of faint white 
dwarfs in binary systems are missing from surveys, and then we calculate the number of white 
dwarfs expected given the observed number of M dwarfs and a mildly decreasing star formation 
history; these two numbers are in good agreement. Finally, we re-assess the distribution of gas 
in the solar neighborhood. We estimate the amount of HI that is missed due to optical depth in 
the 21 cm line from two independent measurements, and we are able to significantly reduce the 
uncertainty in the column density of HI. We also give improved results for the vertical distributions 
of the gas and stars. 

Beginning with Kapteyn (1922), astronomers have attempted to infer the density of dark 
matter in the solar vicinity, pdM; by analyzing star counts and stellar kinematics. In the modern 
era, this study was rejuvenated by Kuijken k. Gilmore (1989a, 1989b, 1991), who analyzed a 
sample of K dwarfs towards the south Galactic pole. Since then, a variety of techniques have been 
developed and much larger samples of stars have been analyzed, as discussed in the comprehensive 
review by Read (2014). In many cases, analysis of stellar kinematics provides a value of the vertical 
gravitational acceleration, K^. For a disk of constant surface density with a flat rotation curve, 
this is equivalent to a measurement of the total surface density, S, since Kz = — 27rGS. The 
surface density of the Galactic disk varies, however, so this relation does not hold exactly (Kuijken 
k Gilmore 1989a), and we evaluate the deviation. As Read (2014) points out, one obtains a 
better estimate of pdm if one has a prior direct determination of the vertical density distribution 
of baryons, Pb{z). Our determination of the baryonic density distribution can therefore tighten the 
constraints on pdm- As a first step in that direction, we take some existing measurements of the 
total surface density within 1 or 1.1 kpc of the plane and use our result for the baryonic surface 
density, S;,, to infer revised values for poM- 

This paper is organized as follows. We focus in Section 2 on the observed surface density of 
M dwarfs, which are the single largest contributor to the stellar surface density. Section 3 presents 
our determination of the surface density of white dwarfs. In Section 4 we summarize the surface 
densities of stars and stellar remnants, while Section 5 discusses the local interstellar gas. We bring 
this all together in Section 6, where we present our results for the local baryon budget. In Section 7, 
we apply these results to determining the local density of dark matter, and discuss the implications 
for disk dark matter, including the possible existence of a thin disk of dark matter in the Galaxy. 
Our results are summarized in Section 8. Appendix A discusses the variation of metallicity with 
height and shows that this has only a minor effect on the inferred scale heights. In Appendix B we 
discuss the stellar halo and show that halo stars within about 3 kpc of the plane are included in 
our analysis. Finally, in Appendix [C] we determine the relation between the vertical acceleration 
due to gravity and the surface density of matter. 
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2. The Surface Density of M Dwarfs 
2.1. Counting the M Dwarfs 

M dwarfs are the largest single contributor to the stellar mass, so we begin by estimating the 
surface density of such stars in the solar neighborhood. We denote stellar masses, measured in 
solar masses, by m. The mass range of M dwarfs extends from the maximum mass of a brown 
dwarf, m-BD = 0.075 (Burrows et al. 2001), which is almost independent of metallicity, to the 
minimum mass of a K dwarf, ttim-k- Following Gould et al (1996) and Zheng et al. (2001), we 
shall identify M dwarfs as those with absolute visual magnitudes My > 8. For the mass-luminosity 
relation adopted by Reid et al (2002, hereafter RGH02), this corresponds to niM-K = 0.69; for the 
mass-luminosity relation of Kroupa et al. (1993), this corresponds to thm-k = 0.65. We shall take 
the average of these and adopt tum-k = 0.67. This value applies in the solar neighborhood; the 
metallicity decreases away from the Galactic plane, and so does the value of niM-K (Appendix!^. 

Gould et al. (1996) used HST observations of distant field M dwarfs to make the first deter¬ 
mination of the surface density of M dwarfs in the solar neighborhood. They inferred My from 
a color-magnitude relation that was independent of height. Binaries were not resolved in their 
observations, so they determined “effective masses”, that is, stellar masses derived by assuming 
that binaries are single stars with the luminosity of the binary. They showed that the stellar mass 
surface densities inferred from the “effective mass function” agree reasonably well with the true 
surface densities. They used the Henry & McGarthy (1993) mass-luminosity relation, for which the 
maximum mass of an M dwarf (i.e., the dividing line between M and K stars) is tum-k = 0.66, 
which is very close to the value we have adopted for local stars, tum-k = 0.67. They found that 
simple one-component models for the number density of M dwarfs, as a function of distance 

above the Galactic plane, z, such as oc secli^ (z/h) (expected for an isothermal, self-gravitating 
disk) or oc exp{—z/h) (expected for isothermal stars in a constant gravitational field), are in¬ 
consistent with the observations. (Throughout this paper we take ^ > 0 and assume that the matter 
distribution in the Galaxy is symmetric about the Galactic plane.) However, simple two-component 
models, such as a double exponential or the sum of a sech^ term and an exponential, 

z 

~h2 

are consistent with the data. Here n*^M,o is the density of M dwarfs in the midplane, (3 is the 
fraction of the midplane density associated with the thick disk, and /12 is taken to be greater than 
hi- They emphasized that the two terms do not necessarily correspond to two physically distinct 
populations, and indeed Bovy, Rix, & Hogg (2012b) have shown that there is no distinct thick disk; 
instead, populations defined on the basis of metallicity have a smooth distribution of scale heights. 
Nonetheless, for simplicity we shall refer to these two terms as describing the thin disk and thick 
disk, respectively. The effective scale height for this number-density distribution (or for the sum of 





(1 - /3)sech ( — ) -I- /3exp 
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two exponentials) is 

/iM = X^^ = (l-/3)/il+/3^2, (2) 

2n*_M,0 

where is the number of M dwarfs per unit area in the Galactic disk. Under their assumption 
that there is no metallicity gradient with height, the scale height hyi is the same as the mass scale 
height, 

/im,M = 5 ]*,m/(2/?*,M,o)- (3) 

For a distribution of the form in Equation ([T]), they inferred that the local surface density of M 
dwarfs is = 12.4± 1.9 Mq pc“^, corresponding to an effective scale height hm,M = 390 pc. For 
the double-exponential fit, they inferred a somewhat larger surface density, = 14.5 ± 2.4 Mq 
pc“^, but a considerably smaller effective scale height, hm,M = 236 pc. 

As pointed out by Kuijken k. Gilmore (1989b), Gould et al. (1996), and ZOl, stellar surveys are 
strongly biased towards stars at large distances from the Galactic plane. Since the functional form 
of the density distribution in the disk is unknown, this leads to a fundamental problem: Even with 
an accurate determination of the stellar density away from the plane, it is difficult to accurately 
infer the density near the midplane from such surveys, and correspondingly to accurately measure 
the stellar surface density. Gould et al. (1996) used a clever trick to mitigate this problem: They 
fit the data to both a sech^ -t- exponential (Equation ([1])) and to a double exponential. Since these 
fits have different midplane densities, they were able to linearly combine them so as to match the 
observed midplane density. It is important to note that although there are significant uncertainties 
in the distances to the stars, and therefore in the stellar density, the uncertainty in the number of 
observed stars is only statistical and is small. One can show that the number of observed stars is 
the same for the two fits for the density, and the linear combination preserves this number. Gould 
et al. (1996) inferred that the contribution of the double exponential to the total density was 
negligible by comparing their results with the number density of early M dwarfs (8.5 < My < 12.5) 
observed locally by Wielen, Jahreiss, k Kruger (1983)-i.e., the sech^-exponential fit is in much 
better agreement with the data than the double exponential fit. These results were revised slightly 
by Gould et al. (1997). 

ZOl extended this work with more data. They considered two cases: one with a constant 
color-magnitude relation, GMR(I), corresponding to no variation of metallicity with distance from 
the plane, and one with a z-dependent relation, GMR(II), which is more realistic (e.g., Bochanski 
et al. 2010; we use roman numerals to distinguish these CMRs in order to avoid confusion with 
the notation for the thin and thick disks). For GMR(II), they also included an estimate of the 
z-dependence of the mean mass-luminosity relation; a decrease in the metallicity results in a lower 
mass for a given luminosity. In order to eliminate contamination by spheroid stars, they considered 
only stars within a distance Zmax of the plane, with Zmax = 3200 pc for GMR(I) and Zmax = 
2400 pc for GMR(II). For GMR(I), they found = 14.3 ± 1.3 Mq pc“^ for the value of = 
0.0180Mq pc“^ inferred from Wielen et al. (1983), corresponding to hm,M = 397 pc. For the 
z-dependent case, they used the ^-dependent GMR(II) to infer the luminosities and distances of 
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the stars. They adopted a z-dependent mass-luminosity relation, but rather than using that for 
each star, they applied it with the mean observed height in each magnitude bin. The lower masses 
associated with lower metallicity led to a reduced midplane density, p*,m,o = O.O153M0 pc In 
this case they found = 12.2±1.6Mq pc“^, which corresponds to an effective mass scale height 
hm,M = 399 pc, almost identical to that for CMR(I). 

Subsequently, the properties of the stellar disk have been inferred from SDSS data. Juric et 
al. (2008) inferred the properties of the stellar disk from observations of M1-M3 stars in a specified 
range of r — i colors, corresponding to a range of absolute magnitudes that was independent of 
distance from the plane. Using a double-exponential fit, they found hi = 300 ± 60 pc, /i 2 = 
900 ± 180 pc, and jS = 0.107, corresponding to /im = 364 pc. In an independent analysis of 
SDSS data, Bochanski et al (2010) also used a double-exponential fit and found hi = 300 ± 15 pc, 
/i 2 = 2100 ± 700 pc, and /3 = 0.04. Even though these properties of the thick disk are very 
different from those found by Juric et al. (2008), the effective scale height is very nearly the same: 
/iM = 372 pc. As noted by Bovy &: Rix (2013), the mass-weighted scale height of all the mono¬ 
abundance populations of Bovy et al. (2012b,c) is 400 pc. Furthermore, Bovy &:Rix (2013) inferred 
a scale height of 370 ± 60 pc All these values are in good agreement with the results of Gould et 
al. (1996) and ZOl. 


2.2. Revising the Results of Zheng et al. (2001) 

The measurement of the surface density of M dwarfs by ZOl can be improved since they used 
pre-Hipparcos data from Wielen et al. (1983). Jahreiss & Wielen (1997) and RGH02 determined 
the local stellar luminosity function with the aid of Hipparcos data, although it should be noted 
that most of the data on M dwarfs still comes from ground-based dataQ Using the mass-luminosity 
relation proposed by RGH02, which is based on the results of Delfosse et al. (2000) for My > 10 
and of Andersen (1991) for My < 10, we find that the local mass density of M dwarfs is 

/3*,M,o = 0.0216 Mq pc“^ (4) 

from the RGH02 data (the earlier data of Jarhreiss & Wielen (1997) gives 0.0198Mq pc“^). This 
density includes the contribution of the M-dwarf companions of stars more massive than M dwarfs, 
which ZOl did not include; the data of RGH02 imply that this contribution is O.OOI 2 M 0 pc“^, 
so that the density of M dwarfs excluding this contribution is about O.O 2 O 4 M 0 pc“^, somewhat 
greater than the value of O.OI 8 OM 0 pc“^ adopted by ZOl for GMR(I). 

We could follow a procedure similar to that of Gould et al. (1996) and ZOl to determine 


^Note that the more recent data imply that the mass density of stars with 8.5 < Mv < 12.5, which ZOl used 
to normalize to the observations, is smaller than found by Wielen et al. (1983); however, the total mass density of 
M dwarfs and their companions (i.e., excluding M-dwarf companions of more massive stars) adopted by ZOl from 
Wielen et al. (1983) is very nearly the same as that from Jahreiss & Wielen (1997). 
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the M-dwarf surface density: combine the sech^-exponential and double-exponential fits so as to 
agree with the local density of stellar systems observed by RGH02, and then use the mean mass per 
stellar system. However, the much more straightforward approach of simply using the effective scale 
height determined by ZOl and replacing the local density that they used with the value measured by 
RGH02 gives very nearly the same answer, so we shall use that. For both GMRs, ZOl found a scale 
height /iM — 400 pc. This scale height is somewhat larger than the M-dwarf scale heights found by 
Juric et al. (2008) and Bochanski et al. (2010), but this is consistent with the finding by Gould 
et al. (1996) and by ZOl that a sech^-exponential fit gives a larger effective scale height than the 
double exponential fit used by the former authors. A virtue of the fit we have adopted from ZOl is 
that it reduces the magnitude of the density gradient at the midplane, which is expected to vanish: 
For example, for the ZOl CMR(II) normalized density ht (i.e., the ht combining a sech^-exponential 
and a double exponential), the relative gradient at the midplane, {d/dz){n^^y[/n^,^M,o)-, is 5 times 
smaller than for the Juric et al. (2008) or Bochanski et al. (2010) double exponentials. 

In Appendix Owe show that the mass-weighted scale height, that is needed to infer the 

surface density is very close to hyi- Given the uncertainties, we shall adopt hm,M — — 400 pc 

based on the results of ZOl. With this value of the effective scale height and the local mass density 
of M dwarfs determined by RGH02, we conclude that the surface density of M dwarfs is 

Z;*,M = = 17.3 ± 2.3 Mq pc“^, (5) 

where we have estimated the error based on the ZOl result for the metallicity-dependent CMR(II). 
This value of the surface density is significantly higher than the value 14.3 Mq pc“^ of ZOl with 
CMR(I) and higher still than the value 12.2 Mq pc“^ of ZOl with GMR(II) because it has a higher 
midplane M-dwarf density and it includes the O.OMq pc“^ associated with the M-dwarf companions 
of more massive stars. 

Our inferred value for the surface density of M dwarfs is the same for the two GMRs since 
we take the midplane density from observation and since the scale heights obtained by ZOl are 
almost exactly the same for the two GMRs. Why is it that ZOl found that the metallicity gradient 
significantly reduced the inferred surface density, from 14.3 Mq pc“^ to 12.2 Mq pc“^? The primary 
reason for this is that they used the average observed height of the stars to infer the effect of the 
metallicity gradient, whereas we use the average aetual height. As they point out, their observations 
are biased toward stars far from the plane; indeed, Figure 2 in ZOl shows that the average observed 
height of the stars in their sample is about 1100 pc. On the other hand, the actual value of the 
average height of the stars is of order the scale height, 400 pc, which is much less, thereby reducing 
the effect of the metallicity gradient. A second difference between our analysis and theirs is that 
they found that the metallicity variation had a substantial effect on the mass function, whereas 
our assumption of a universal IMF means that there is no effect. As discussed in Appendix in 
our work the effect of the metallicity gradient enters only in its reduction of the maximum mass 
of an M dwarf niM-Kiz) with height. The assumption that the IMF is independent of metallicity 
is supported by the analysis of the mass functions of globular clusters, which have a wide range of 
metallicity but have mass functions that depend primarily on the ratio of the time for the cluster to 


dissolve by two-body encounters to the age of the cluster (de Marchi, Paresce, &: Portegies Zwart 
2010 ), not on the metallicity. 

We are not aware of other direct determinations of the M-dwarf surface density that we can 
compare to our result. We shall compare our result for the total stellar surface density with other 
results in Section [4.41 


3. White dwarfs 

White dwarfs make a significant, but uncertain, contribution to the surface density. FOG 
combined white dwarfs, neutron stars and black holes into a single category, but as we shall see 
below, this category is completely dominated by the white dwarfs. We infer that they found a total 
surface density for these objects of 6.9 Mq pc“^, which is is larger than the values given in the 
references below. For white dwarfs, Sion et al. (2009) find a local density of 4.9 x 10“^ WD pc“^, 
which is in good agreement with the previous estimate of (5 ± 0.7) x 10~^ WD pc“^ by Holberg et 
al. (2002), although it is somewhat less than the estimate of 6.0 x 10“^ WD pc“^ by Reid (2005). 
Assuming a scale height hi = 325 pc, Sion et al (2009) estimated that the surface density of white 
dwarfs is W,WD = 3.2 WD pc“^. They argued that there is no thick-disk contribution for white 
dwarfs, but this is difficult to accept in view of the well-established thick disk for the MS stars that 
are the precursors of white dwarfs (Juric et al. 2008, Bochanski et al. 2010). We note also that 
the vertical velocity dispersion of the DC and DQ white dwarfs in their sample exceeds 37 km s“^, 
which corresponds to a scale height significantly greater than 325 pc. Based on observations that 
extend well beyond the 20 pc local sample, a number of authors have presented direct evidence for 
a thick-disk component, with (5 = 0.2 (Reid 2005), /? = 0.16 (Rowell & Hambly 2011), and /3 = 0.33 
(Fuhrmann et al. 2012). If we combine the anomalously large halo component found by Rowell 
& Hambly with their thick-disk component, their value of /3 would become 0.21. It is reasonable 
that (3 is larger for white dwarfs than it is for main sequence stars since on average the precursors 
to white dwarfs are older than main sequence stars. Since the white dwarfs observed in all these 
cases are almost all within 300 pc, it is not possible to infer the scale heights directly from the 
observations. The scale heights determined by Juric et al. (2008) give h = 420 pc for /3 = 0.2 and 
h = 500 pc for j3 = 0.33. In Section [321 we infer that the scale height of white dwarfs is about 430 
pc, close to the Juric et al. value for (3 = 0.2. For the Sion et al. (2009) density, an effective scale 
height of 430 pc gives a surface density by number of W,WD = 4.2 WD pc“^. 

Recently, Katz et al. (2014) have suggested that the number of white dwarfs could be under¬ 
estimated by up to a factor 2 based on the difficulty in detecting white dwarfs in binary systems. 
They point out that the number of bright white dwarfs with main sequence companions is about 
equal to the number of bright single white dwarfs, but the number of faint white dwarfs with such 
companions is much less than the number of faint singles. They argue that the number of white 
dwarfs in binaries is likely to be at least as large as the number of single white dwarfs, which leads 
to their estimate that the true number of white dwarfs has been underestimated by about a factor 
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of two. A simple way of expressing their result is that if n*^wD,s is the density of single white dwarfs 
and /wD,s is the fraction of white dwarfs that are single, then the total density of white dwarfs is 


n*,wD 


n*,wD,s 

fwD,s 


( 6 ) 


Sion et al. (2009) find 89 isolated white dwarfs within 20 pc and argue that this is 80% complete, 
which implies that locally the density of single white dwarfs is n* ,WD,s = 3.3 X 10 ^ WD pc 
Katz et al. (2014) infer /wD,s = 0.4 from Raghavan et al. (2010), which implies that the total 
density of white dwarfs is n*,wD = 8.3 x 10“^ WD pc“^, which is 1.7 times larger than the value 
of 4.9 X 10“^ WD pc“^ cited by Sion et al. (2009). For a white-dwarf scale height of 430 pc, this 
corresponds to a surface density of A'^^wd = 7.1 WD pc“^. 


In view of the substantial discrepancy between the white-dwarf density found by Sion et al 
(2009) and that we infer from Katz et al. (2014), we now provide an alternative estimate of the 
white-dwarf density that is based on the IMF and is independent of binarity. 


3.1. Ratio of White Dwarfs to M dwarfs 

We begin by calculating the surface density of white dwarfs produced by the stars that have 
formed in the disk. We describe the IMF by ipim), where ^lJ{m)dlnm is the fraction of stars 
(including brown dwarfs) born between m and m + dm and where m is measured in solar masses. 
The rate at which stars are born per unit area in the mass range dm is then 

dN^{m) = N^,T'4’{m)d\nm, (7) 

where is the total rate at which stars are born per unit area. The star formation rate (SFR) 
per unit area in solar masses per unit time is S* = A*'7’(m), where (m) ~ 0.5 is the average mass 
in the IMF. 

Most IMFs (e.g., Kroupa (2002), Chabrier (2005), and PMHlljl have a power-law form, 
'll: oc m~^, for mu > m > 1, where F = 1.7 for the Kroupa (2002) individual star IMF and F = 1.35 
for the other two IMFs; m^, the upper bound on the stellar mass in the power-law regime of the 
IMF, exceeds 100 Mq. In this power-law regime, the stellar birthrate per unit area can be written 
as 

dN^,{m) = RS*m~'"(ilnm, (8) 

where B is a constant that depends on the IMF. One way to characterize IMFs is through the 
total mass of stars formed per high-mass star, We define a “high-mass star” as one with 


^The Kroupa (2002) IMF is a broken power law, '0 oc with 4 segments; the Chabrier (2005) IMF is a log 
normal function of m for m < 1 and a power law for m > 1; and the PMHll IMF has the form oc m~^[l — 
exp —(m/mch)^^'"], which smoothly transitions from ip oc m'*' at low mass to ip oc m~^ at high mass. We have 
changed notation from PMHll, replacing c«(m) with dNt,{'m)/dlnm. 
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niu > m > rrih, and choose nih = 8 so that the number of high-mass stars corresponds to the 
number of core-collapse supernovae, neglecting the complications associated with binary evolution 
(Sana et al. 2011). The rate of high-mass star formation per unit area is 


N*h — 


^dlnm = 




Tml 


where 


</>/), = 1 - 


mh\ 

mu) 


is close to unity. The mass of stars formed per high-mass star is so that 


B = 


d'h4>h 


(9) 


( 10 ) 


( 11 ) 


The quantity B is relatively constant for the three IMFs we are considering: B = (0.276, 0.255, 0.237) 
for the Kroupa (2002), Chabrier (2005) and PMHll individual star IMFs, respectively, whereas Hh 
varies by more than a factor 2; fih = (213, 90, 97). 

We describe the star-formation history in terms of the ratio of the star-formation rate to the 
value averaged over the age of the disk at the solar circle, to (cf. Miller & Scalo 1979), 


b{t) = 


{N*t) ’ 


where 


(77*t) — — f ^*T{t) dt. 
to Jo 


( 12 ) 


(13) 


Note that to corresponds to the present time; it is the time that has elapsed since the onset of star 
formation in the solar vicinity. We further note that by definition we have 

rio 

(14) 


1 

— / b{t) dt = 1. 

-0 Jo 


to Jo 

The birthrate per unit area of white dwarfs at time t due to stars which had a main-sequence 

is 


mass in the mass range rupre to rupre + dw-pre 

diVwD(m.pre,i) = WT(tpre)V’("ipre)dlnmpi.e, 


(15) 


where tpre = t — Tpre IS the time at which the precursor was born and Tpre is the lifetime of the 
white-dwarf precursor of mass rupre, including post main-sequence evolution. The lowest-mass 
white dwarf precursor has a precursor lifetime equal to the age of the disk, to', label this mass 
"Jpre,o- Under the assumption that white dwarf precursors have a maximum mass of 8Mq, the 
total surface density of white dwarfs is then 

rS ^pre 

/ ?Tipj.g d In TTlpre / f^(tpre)dtpre) ( 16 ) 

= aiB{t^)to, (17) 
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where (S*) is the star-formation rate averaged over the time to and ai is a dimensionless parameter 
of order mpre,o~'") which as we shall see is of order unity. Adopting a power-law form for the 
precursor lifetime as a function of mass, 


"^pre — 


/ ?7ipre,0 

V ^pre 


i 

to, 


(18) 


we find that in the mass range 0.9 < rupre < 1.6, setting i = 3.5 and mpre,o = (1-06, 0.99) for 
metallicities of Z = (0.019, 0.008), respectively, gives an accuracy better than ~ 6% based on the 
evolutionary tracks of Girardi et al. (2000). Given the weak dependence of ?n,pre,o on metallicity, 
we shall neglect it and set mpre,o = 1 in our numerical evaluations. 

As a simple model for a variation in the star formation rate, we adopt a linear form for b{t) 
(Eq. m, 

b{t) = 2 - 6o - 2(1 - bo)t/to, (19) 

where bo = b{to) is the current value of b{t). Then the coefficient ai in Equation (fT7|) can be 
evaluated as 

~ mt + 2T-rb„) _^ 

r(r + f)(r + 2<)mp,.y rgr’ ' ' 

where we have ignored terms of order 8“^, and we have adopted £ ~ 3.5 for 1.6 < rupre < 8 since 
that mass range does not contribute much to the integral. Eor a constant star formation rate 
{bo = 1) and mpi.e,o — this gives oi = 0.49, whereas in the opposite case in which the current 
SFR is very small {bo ~ 0) we have ai = 0.58. 


We now compare the predicted surface density of white dwarfs with the surface density of M 
dwarfs, 

E*,M = (21) 


where m is the M-star mass fraction produced by a given IMF. In writing this equation, we have 
ignored the weak dependence of Fm, m on metallicity; the effect of the variation of mean metallicity 
with height is discussed in Appendix^ For the Kroupa (2002), Ghabrier (2005) and PMHll IMFs, 
Fm,M is given by Fm ,m = (0.45, 0.29, 0.33), respectively. Inserting Equation (f2T]l into Equation 
(ffB, we find that the ratio of the surface density of white dwarfs by number to the surface density 
of M dwarfs by mass is 

^ ( 22 ) 

Alternatively, in terms of the midplane densities, we have 

?T-WD,0 _ / QiF 
P*,M,0 \Fm, M 


hm,M 

hwD 


(23) 


Although the white dwarf scale height, /iwo, is not directly measured, its value relative to that for 
M stars can be inferred, as we shall see below. 










- 12 - 


3.2. The Mean Age and Scale Height of White Dwarfs 


It remains to estimate the ratio of the scale heights of white dwarfs and M dwarfs. The 
velocity dispersion of stars is observed to increase with age; from an analysis of the 16682 nearby 
F and G dwarfs in the Geneva-Copenhagen Survey of the Solar Neighborhood (Nordstrom et 
al 2004), Holmberg et al. (2009) found that the vertical velocity dispersion, aw, increases as 
g^ggO.53 _ _ ^^'jO.53^ where G is the time at which the star was born, so that to — t* is its age. We 

find that the results of F06 imply that the scale height, h, scales as for stars with My > 3, 
so that h (X {to — For a white dwarf, the relevant age is that of its precursor, to tp 

Averaging over all the white dwarfs and M dwarfs, we have 

(t0-tpre)"°-"' 


‘^pre- 


^WD 

hu 


(^0 ^*,m) 


(24) 


For the form for b{t) that drops linearly from t = 0 to t = to that we have adopted (Eq. [19]), 
the average birth time of M dwarfs is 


(G,m) = 


1 


1 + 2^0 ) to, 


(25) 


which corresponds to an average age of ^to for a constant SFR {bo = 1) and |to for 6o = 0 (recall 
bo = b{to))- The average birth time of the stars that became white dwarfs is 

^^0 "^pre 


{tpve) — 


1 


Awd 
B{t,{t)) 
A'wd 
a2R(il*(t))t§ 


pre J ‘'predtpj^g, 


f f 

/ dlnm 

J ^^nre,0 " 0 

/■8 dlnmpre 

/ r 7 ?r / Ktpre)t 

-'?Tlpre,0 ''^pre -'0 


N^{m, 

^pre)^; 

^0 "^pre 


prey^^pre'-^f'pre? 


dtr. 


A'wd 


where 


02^0 

ai 


e 


0-2 


2 + 


{i - r)bo 

r + £ 


i+ii« 


1 


(26) 

(27) 

(28) 

(29) 

(30) 


F(F + 2^)(r + 3£)mpre,o^ [ T + i \ 3 V 2 7 FS^ ’ 

where again we have ignored terms of order 8“^ and have adopted £ ~ 3.5 for 1.6 < rupre < 8. 
For a constant SFR and mpre,o — Ij this corresponds to an average age of white-dwarf precursors 
of 0.59to, whereas for bo = 0, the average age is O.Tlto- The ratio of the age of the white-dwarf 
precursors to that of the M dwarfs is 


{to tpre) 3 Oi 02 

{to-t*,M) 2 oi ( 1 - 46 o )’ 


(31) 


which is 1.18 for bo = 1 and 1.06 for bo = 0—i.e., 1.12 ± 0.06 for 1 > 6o > 0. Normalizing to the 
case bo = ^, which is consistent with the results of Aumer & Binney (2009) and approximately the 
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value we get in Parravano, McKee & Hollenbach (in preparation), we find that Equations (j24l) and 
(I3ip imply that for P = 1.35 the ratio of the scale heights is 

1.085 -0.080 (32) 

um \2 J 

to within 0.1%. In Section YI?2[ we adopt hyi = 400 pc for the M-dwarf scale height; for bo = 
this gives /iwD = 434 pc. 


3.3. Surface and Volume Densities of White Dwarfs 


In terms of the observed surface density of M dwarfs, the predicted surface density of white 
dwarfs for the PMHll IMF is (Eq. 


VwD = 


6.7 +1.1 ( i-6o 




17.3 Mq pc 


-2 


WD pc 


-2 


(33) 


The local volume density of white dwarfs is (Eq. 


?T'WD,0 = 


7.6 X 10"^ + 1.8 X 10“^ ( ^ - 6r 


P*,M,0 


0.0216 Mq pc 


-3 


WD pc-+ 


(34) 


For the Chabrier (2005) individual-star IMF, V-wd and n-wD are increased by a factor 1.21. 


We conclude that if the star formation rate has decreased linearly in time over the life of the 
disk so that it is now half the average rate {bo = 0.5), then the PMHll IMF implies that the density 
of white dwarfs is u,*^wD = 7.6 x 10“^ pc~^, and the Chabrier (2005) individual-star IMF implies 
that it is n*^wD = 9-3 x 10“^ pc“^. Equation (l3H) shows that this density estimate is not sensitive 
to the exact value of the current rate divided by the average rate. We shall adopt the average of 
these two values, n*,wD = 8.5 x 10“^ pc~^, which is very close to the value of 8.3 x 10“^ WD 
pc“^ that we inferred from Katz et al. (2014). This analysis confirms their conclusion that the 
number density of white dwarfs has been significantly underestimated by a number of workers. To 
obtain a rough estimate of the error in these numbers we add in quadrature the uncertainty in the 
column density of M dwarfs (2.3/17.3) and the uncertainty in the IMF, which we take to be 10% 
based on the difference between our adopted density and the values implied by either the PMHll 
or Chabrier IMFs. The result is an uncertainty of 17%. 


To convert these results for number densities to mass densities, we adopt the Holberg et al. 
(2008) value for the mean mass of local white dwarfs, 0.665 Mq. This value is inferred from a 
catalog of white dwarfs within 20 pc of the Sun, which they infer to be 80% complete; this catalog 
is almost identical to the later catalog of Sion et al. (2009). The masses are determined from the 
stellar surface gravities and effective temperatures; in addition, most of the stars have trigonometric 
parallaxes. The masses of the white dwarfs with well determined masses (uncertainties less than 
10%) in this catalog range from 0.40 Mq to 1.25 Mq. Over half the stars have mass uncertainties 
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less than 0.04 Mq. With this value of the mean mass of white dwarfs, we find that the local mass 
volume density is /0*o,WD = 0.0056 ± 0.0010 Mq pc“^ and that the local mass surface density is 
S*,WD = 4.9 =h 0.8 Mq pc~^. 


4. Stellar Surface Densities 

The standard reference for stellar surface densities in the solar neighborhood is F06, which 
builds on the previous work by Holmberg & Flynn (2000, 2004). These authors assumed that 
each stellar component was isothermal and related the local mass density to the surface density 
by solving the Poisson-Boltzmann equations. For M dwarfs, they took the surface density, 
from the work of Gould et ah (1998) and then inferred the local mass density of these stars, 
whereas for all other stars they inferred the surface density from the observed local density 
and the calculated scale height. Holmberg &: Flynn (2004) introduced a thick stellar disk in their 
model for the vertical distribution of stars. FOG updated this model, and from their results we 
infer that the thick disk has 10% of the local volume density of stars (/3 = 0.1) and a scale height 
/i 2 = 1000 pc; they assumed that main-sequence stars with My < 4 were too young to have a thick- 
disk component. FOG find 28.5 Mq pc“^ in the thin disk and 7 Mq pc“^ in the thick disk. They 
do not specify the composition of the thick disk, although it is clear that it has no main sequence 
stars with My < 4. We assume that the other stars in the thick disk, which are all long-lived, have 
the same proportions as in the thin disk. (Since the stars in the thick disk are older than those in 
the thin disk, the proportion of red giants and white dwarfs in the thick disk should be somewhat 
larger, but we omit that complication here.) After distributing long-lived stars from the thick disk 
in proportion to their abundance in the thin disk, we represent the FOG model as shown in Table 
[TJ FOG include halo stars as a separate category. As discussed in Appendix iBl determinations of 
the surface density of stars generally consider stars below ~ 3 kpc; if halo stars are not included as 
a separate category, they are included in the thick disk. About 2/3 of the stars in FOG’s halo are 
within 3 kpc of the plane, and they are included in the table. The value of the local surface density 
of stars and stellar objects found by FOG is 36.1 Mq pc“^. It is of interest to compare their result 
for main sequence stars to that found many years ago by Miller and Scalo (1979). These authors 
found a surface density of 27 Mq pc“^; since they estimated the scale height of these stars as 325 
pc, this does not include halo stars. The FOG value for main-sequence stars in the disk (in Table [H 
visible stars minus giants and halo stars), 25.8 Mq pc“^, is quite close to this. As we shall see, our 
value for this quantity, 26.4 Mq pc“^, is even closer to the Miller & Scalo (1979) value. 


4.1. Main-sequence stars and giants 

We update the FOG results in several ways. First, we use the local densities of main-sequence 
stars measured by RGH02 (see Table [T]) . The local density of main sequence stars we calculate 
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Table 1. Local Stellar Surface Densities (Excluding Halo Stars with z > 3 kpc)'^ 


Description 

S*(F0G)^ 

P*o 

h 


F*i.i 


{Mq pC-2) 

{Mq pC-3) 

(pc) 

(Mq pC“2) 

{Mq pC“2) 


My < 3 

1.5 

0.0018 

140 

0.5 

0.5 

3 < My < 4 

1.1 

0.0018 

236 

0.8 

0.8 

4 < My < 5 

2.2 

0.0029 

384 

2.2 

2.1 

5 < My <8 

7.2 

0.0072 

400 

5.8 

5.4 

8 < My (M dwarfs) 

13.8 

0.0216 

400 

17.3 

16.2 

Giants 

0.5 

0.0006 

344 

0.4 

0.4 

Halo ( 2 ; < 3 kpc)'^ 

0.4 



0 

0 

Visible stars 

26.7 

0.036 


27.0 

25.4 

Brown dwarfs (BD) 

2.3 

0.0015 

400 

1.2 

1.1 

White dwarfs (WD) 

6.9^ 

0.0056 

430 

4.9 

4.5 

Neutron stars (NS) 


0.0001 


0.2 

0.1 

Black holes (BH) 


0.0001 

400® 

0.1 

0.1 

Total 

35.9 

0.043 


33.4 ± 3 

31.2 ±3 


is the local stellar surface density based on F06. Numbers to the right of the vertical line 
are based on our analysis: /9*o is the local stellar mass density in the midplane, h is the effective 
scale height, h = S*/(2p*o)) and is the local surface density of stars within 1.1 kpc of the 

midplane. 

^The 7 Mq pc“^ in the FOG thick disk has been distributed to the different categories as described 
in the text. The FOG halo stars above 3 kpc (0.2 Mq pc“^) have been excluded. 

"^Halo stars within about 3 kpc of the Galactic plane are included in our entries for stars and 
stellar remnants (Appendix [B]) . 

^This value includes neutron stars and black holes. 

®The scale height of black holes is unknown and has been assumed equal to that of M dwarfs. 
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from the RGH02 results is /3*o(MS) = 0.035 Mq pc which is the same as that from Foelfl adding 
0.0006 Mq pc“^ for giants then brings the total local density of visible stars up to 0.036 Mq pc“^. 
The local density of M dwarfs measured by RGH02 is larger than that adopted by F06, and as a 
result the M-dwarf surface density we determined in Section [221 = 17.3 Mq pc“^, is larger 

than the value we infer from the F06 model, S* m = 13.8Mq pc The RGH02 density for stars 
with My < 3 is significantly less than that used by F06. However, in their Table 3, F06 cite a 
subsequent, more complete analysis of these stars that is within 10% of the RGH02 value, so we 
regard the value we adopt as more accurate. 

Next, we reduce the F06 surface densities for stars earlier than M dwarfs by adjusting their scale 
heights to be consistent with more recently determined values. For low-mass stars with My > 4, 
F06 found that the scale height of the thin disk was hi ~ 400 pc. However, as discussed above, two 
subsequent analyses have found that in fact hi ~ 300 pc (Juric et al. 2008, Bochanski et al. 2010). 
Other than the larger value of the thin-disk scale height in their model, it is very similar to the 
model of Juric et al. (2008); /? = 0.1 vs. j3 = 0.107 and /i 2 = 1000 pc vs. /12 = 900 pc. Although 
Bochanski et al. (2010) found a very different model for the thick disk than Juric et al. (2008), 
the effective scale heights were almost the same: h = 364 pc for Juric et al. (2008) and h = 372 pc 
for Bochanski et al. (2010). By contrast, the effective scale height for stars with My > 5 implied 
by the F06 model \s h = (0.9 x 405 -|- 0.1 x 1000) pc = 464 pc. Previous authors had also found 
smaller scale heights than F06: Recall that ZOl found h ~ 400 pc; in another case, RGH02 cite 
evidence that the scale height for stars with My < 3 is only 100 pc, whereas the F06 scale height 
for stars with My < 2.5 is 145 pc. A strength of the F06 analysis is that it is based on solving 
the Poisson-Boltzmann equation for each type of star separately. We therefore assume that the 
relative values of the F06 scale heights are correct, and adjust the overall normalization. To be 
consistent with our analysis for M dwarfs, we adopt h = 400 pc for the scale height of stars with 
My > 5 and reduce the effective scale heights for stars with My < 5 by a factor 400/464 = 0.86. 
The resulting scale heights are given in Table [U we have followed F06 in assuming that the thick 
disk has no main-sequence stars with My < 4. To obtain the surface densities, we combine these 
revised scale heights with the stellar densities from RGH02. The corrected scale height for stars 
with My < 3, which is 140 pc, is significantly greater than the value cited by RGH02, so it is 
possible that the surface density of these stars is over-estimated. For giant stars, we use the value 
of the local density given by F06; the reduction in the effective scale height for these stars reduces 
their surface density from 0.5 Mq pc“^ to 0.4 Mq pc“^. Overall, since our value for the M-dwarf 
surface density is larger than the value F06 inferred from Gould et al. (1998), our value for the 
surface density of main-sequence stars, = 26.7 Mq pc“^, is slightly higher than their value, 

5^*,MS = 25.8 Mq pc“^. The F06 model also includes a small contribution from the stellar halo; as 
discussed in Appendix [Bl halo stars within about 3 kpc of the Galactic plane are included in our 
results for stars and stellar remnants in Table [TJ 


®RGH02 quote a lower value, however: p*o = 0.031 Mq pc 
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4.2. Brown dwarfs. 

The density of brown dwarfs is uncertain, since the fraction of stars that are brown dwarfs 
is inferred to be smaller in the field than in young clusters. In PMHll, we determined that 
observations of young clusters imply that the number ratio of low-mass main-sequence stars (0.08 < 
m < 1.0) to stars (mainly brown dwarfs) in the mass range 0.03 < m < 0.08 is = 3.76, where 
the superscript indicates that this number counts the individual stars in binaries. On the other 
hand, the ratio of the number of main-sequence stars with m < 1 (almost all of which have 
m > 0.08) to the number of brown dwarfs in the Kirkpatrick et al. (2012) catalog of stars within 
8 pc (almost all of which are in the mass range 0.03 < m < 0.075) is 5.8, corresponding to fewer 
brown dwarfs. In order to infer the mass density of brown dwarfs, it is also necessary to know the 
distribution of brown-dwarf masses, which is often assumed to be a power law, dM^^/dm oc m““, 
where is the number of stars with masses less than m. In the limit of small masses, the IMF 

of PMHll approaches the form dM^t/d\nm oc m'*', so that a = 1 — 7. In PMHll, we estimated 
7 = 0.51, corresponding to a = 0.49; this is in good agreement with the value a = 0.511q 27 
determined from microlensing observations by Sumi et al. (2011). The power law inferred from the 
microlensing obsevations is also consistent with the Kroupa (2002) power law in this mass range, 
7 = 0.7 ± 0.7. However, Kirkpatrick et al. (2012) infer a ~ —1, corresponding to 7 ~ 2. The 
observational estimates for the density of brown dwarfs are thus uncertain due to the apparent 
contradiction between the values obtained from young clusters and microlensing on the one hand 
and those obtained from infrared observations of the field on the other. We adopt the value 7 = 0.51 
given by the PMHll IMF, which is consistent with observations of brown dwarfs in young clusters, 
with the microlensing observations, and with Kroupa’s (2002) value. For this choice, the total mass 
surface density of brown dwarfs is about S*,bd = 1.2 Mq pc“^ and the mass surface density at 
2; < 1.1 kpc is 1.1 Mq pc“^, but with considerable uncertainty in both values. For example, if we 
ignored the result of Sumi et al. (2011) and adopted the result of Kirkpatrick et al. (2012) for 
the density of brown dwarfs, then these values of the surface density would be reduced by a factor 
3.76/5.8=0.65. We obtained the mass surface density at z < 1.1 kpc by assuming that the scale 
height of BDs is equal to that of M dwarfs, 400 pc, so that 93.6% of the BDs are within 1.1 kpc of 
the mid plane. 


4.3. Neutron Stars and Black Holes 

Neutron stars are often born with high velocities, so it is not possible to directly relate the 
surface density of neutron stars to that of their progenitors. Sartore et al. (2010) have considered a 
range of models for the velocity distribution of neutron stars at birth and have determined the local 
density distribution as a function of distance from the plane. We omit their case IE, which has a 
scale height for neutron stars of only 33 pc. They assumed a fiducial case in which there are a total 
of Nns = 10^ neutron stars in the Galaxy, and their results scale with A^g = ■ Evaluating 

the surface density of neutron stars within 1.1 kpc of the plane in their models and assuming a 
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mean neutron star mass of 1.4Mq, we find S*i,i^ns = (0-3 ± 0.1)Nq Mq pc“^ for the local mass 
surface density of neutron stars within 1.1 kpc of the plane. The average total surface density in the 
Sartore et al. (2010) models we consider is twice this, 5 ]*,ns = (0.6±0.2)iVg Mq pc“^. Ofek (2009) 
performed simulations very similar to Sartore et al. (2010), and the average of his two simulations 
give S*i.i,NS = (0.55 ± 0.05)A^9 Mq pc“^and S*,ns = (0.97 ±0.13)A^9Mq pc“^. Averaging the 
values from Ofek and Sartore et al. (2010), we find S*i.i,ns = (0.4 ± 0.11)A'"9 Mq pc“^and S*,ns = 
(0.8 ± 0.24)A'9 Mq pc“^. Similarly, we average the values of Sartore et al. (2010) and Ofek to 
obtain the local midplane density of neutron stars. We find /9*o,NS = (4.7± 1.6) x 10“^A^9 Mq pc“^. 
Using the PMHll IMF and the Hurley et al (2000) evolution code, which relates the mass of the 
progenitor to the likelihood of producing a neutron star, and scaling to the M-dwarf surface density, 
we estimate Ng = 0.25 (a more detailed description of this calculation will appear in Parravano, 
McKee & Hollenbach, in preparation). 

There are no direct observational estimates of the surface density of black holes, so we infer 
that in an analogous manner to the computation of the number of neutron stars. We use the 
PMHll IMF and the Hurley et al (2000) evolution code, which relates the mass of the progenitor 
to the likelihood and mass of the resultant black hole, and we scale the surface density of black 
holes to the local M dwarf surface density. In order to estimate the volume density of black holes 
in the midplane, we assume that the scale height is similar to that of M dwarfs; the actual value of 
the scale height is unknown. 


4.4. Results and Comparison with Previous Work 

The results on the local stellar surface densities are collected in Table 1. We have estimated 
errors for the M-star surface density (13% from ZOl) and for the surface density of white dwarfs 
(17%, which is dominated by the uncertainty in S*^m)- We lump all the other visible stars together 
and assign an error of 15%; for the brown dwarfs, neutron stars, and black holes, we assign an error 
of 30%. The resulting error in the total surface density is ~ 3 Mq pc“^ and is dominated by the 
uncertainty in the surface density of main sequence stars. Our final value for the local total surface 
density of stars and stellar remnants by mass is then 

S* = 33.4T3 Mq pc-2, (35) 

whereas the local surface density of visible stars (main-sequence stars and giants) is 

Svis = 27.0 ± 2.7 Mq pc-2. (36) 

The surface density of stars and stellar remnants within 1.1 kpc of the Galactic plane is 

S*i.i = 31.2±3 Mq pc-2, (37) 

where we have approximated the error as being the same as that for S*. 
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These results are very close to those of F06: Our results for S* and S*i,i are 8% and 3% 
less than theirs, respectively, whereas our result for Svis is 1% more than theirs. However, the 
excellent agreement between our results and those of F06 is fortuitous, since the mass is distributed 
somewhat differently, but in such a way that the sum of all the components is about the same. We 
find there to be less mass in main-sequence stars earlier than M, in brown dwarfs, and in white 
dwarfs, but more mass in M dwarfs. We have discussed the reasons for these differences in Sections 
12.21 IM dwarfs). 14.11 Imore massive stars), [42] (brown dwarfs), and[3](white dwarfs). This different 
distribution in the mass, especially that of M dwarfs and white dwarfs, is one of the main results 
of this paper. 

Our result for the local surface density of visible stars is also within the errors of the results 
of Bovy et al. (2012b), who found Svis = 30 ± 1 Mq pc“^ for a Chabrier (2001) log-normal IMF. 
Their result has some dependence on the IMF; for example, they find Svis = 29 Mq pc“^ for a 
Chabrier (2003) IMF. They did not take binarity into account in their analysis, but Bovy & Rix 
(2013) show that binarity would not affect their estimates of the surface density of visible stars 
within 1.1 kpc. Furthermore, we note that their result is based on observations of G stars, which 
cover a narrow range of mass. As a result, relatively small errors in the mass-luminosity relation 
can alter the fraction of the stellar mass that is in G stars and thus make a significant contribution 
to the error in the total mass of stars. Their quoted error does not allow for this effect (Bovy, 
private communication). 

Finally, we note that the midplane density of stars is 

p*o = 0.043 ± 0.04 Mq pc-3, (38) 

where we have estimated the error to be 10%; this is considerably larger than the 4% difference 
between our result and that of FOG, but there is a substantial overlap in the data. The effective 
scale height of the stars is then 

K = ^ = 388 ± 40 pc, (39) 

2p*o 

where we have again adopted an error of 10%, which is the difference between the scale height 
found by Juric (2008) et al. and that found by ZOl. 


5. The Local Interstellar Gas 

Having determined the densities of stars and stellar remnants in the solar vicinity, we turn now 
to a discussion of gas in the local interstellar medium-i.e., within a cylindrical radius of 1 kpc of 
the Sun. We consider the molecular (H 2 ), atomic (HI) and ionized (HII) gas in turn. Read (2014) 
and Hessman (2015) have emphasized the uncertainties in existing determinations of the amount 
of local interstellar gas, and our goal in the present section is to reduce those uncertainties. 
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5.1. Molecular Gas: H 2 

In their review, Heyer & Dame (2015) state that the best value for the local surface density 
of H 2 remains that of Dame et al. (1987), who found that the distribution of molecular gas within 
1 kpc could be described by a mean midplane density of H nuclei hH(H 2 ) = 0.2 cm“^ and a half¬ 
width half-maximum of 87 pc. Their density estimate was based on a CO to H 2 conversion factor 
of 2.7 X 10^° cm“^/(K km s“^); using the value of 2 ± 0.6 x 10^*^ cm“^/(K km s“^) suggested by 
Bolatto, Wolfire, h Leroy (2013) reduces the central density to 0.15 cm“^. (It should be noted that 
this empirical conversion factor includes the H 2 that lies in surface layers that have little CO but 
mostly C"'', the so-called “dark” or “hidden” H 2 .) Dame et al. (1987) pointed out that it was not 
possible to determine the shape of the spatial distribution of the local molecular gas because the 
number of massive molecular clouds, which contain most of the mass, is too small. They adopted 
a Gaussian distribution oc exp—( 2 /IO 5 pc)^ based on observations of more distant molecular gas 
by others. The resulting properties for the local molecular gas are given in Tabled In converting 
from the column density of hydrogen to the surface density of mass, we have assumed that the mass 
in helium and heavier elements is 40% of that in hydrogen, so that the mass per hydrogen atom is 
1.4 X mn = 2.34 x 10“^^ g (we shall henceforth refer to this as the He correction). All of this is 
close to the midplane, at 2 ; < 1.1 kpc. The uncertainty is dominated by the CO-to-H 2 conversion, 
which is a factor 1.3. 


5.2. Atomic Gas: HI 

Based on the Hat Creek Survey of the 21 cm emission in the Galaxy (Heiles & Habing 1974), 
Heiles (1976) determined the azimuthally-averaged column density of neutral atomic hydrogen at 
Galactic latitude b for |6| > 10° and for declinations 5 > —30° over the velocity range —92 km s“^< 
u < 75 km s“^. Summing his results for gas above and below the Galactic plane, we find A"^*)™(|6|) = 
[(7.45 ± 0.06) CSC |5| — (2.39 ± 0.14)] x 10^*^ cm“^, where the superscript “thin” indicates that the 
gas is assumed to be optically thin. These observations were contaminated by stray radiation-i.e., 
radiation from directions other than that in which the telescope was pointed. Heiles et al. (1981) 
addressed this problem by reobserving some areas of the sky with the horn reflector antenna at Bell 
Telephone Laboratories, which has much less stray radiation than the Hat Creek telescope. Based 
on their results, we infer that the column density of optically thin HI summed over latitudes +b 
and —6 is 

^ [(-7 45 ^ Q Qg^ 1^1 _ ^3 72 ± 0.14)] X 10^° cm-2. (40) 

We shall refer to this as the “Heiles model” for short. Heiles et al. (1981) interpreted the second 
term as being due to the fact that the Sun is in a low-density region created by stellar winds and 
supernovae. Indeed, the Sun is located in the Local Bubble, which has a density that is much less 
than the average interstellar value (Cox &: Reynolds 1987). We can then write 

A^f"(161)=iV*ji;“cscl61-iVfB^ 


(41) 
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where is the column density of optically thin HI through the entire plane of the Galaxy 

outside the Local Bubble and is the average deficit of optically thin HI across the diameter of 

the Local Bubble. In reality, the Local Bubble is very aspherical, so after averaging over Galactic 
longitude, Afgg™ should depend on Galactic latitude; this dependence is ignored in the Heiles 
model. The average radius of the Local Bubble is i?LB = A^lb°/(^^o^™) ~ 60/nQ'^“ pc, where 
is the midplane density of optically thin HI. The total mass surface density of gas, including He, 
corresponding to Aig*^^ = 7.45 x 10^° cm“^ is = 8.34 Mq pc“^. 

We can check the accuracy of the Heiles model for the HI in the solar neighborhood based on 
the results of Lockman & Gehman (1991; hereafer LG91) from the Bell Laboratories HI Survey 
(Stark et al. 1992). LG91 provide values of Aig')“ weighted by sin|6| for |6| > (45°, 80°) and 
S > —40°. Our comparison is approximate: the Bell Laboratories HI Survey covers a broader range 
of velocities than the Hat Greek Survey, and in our comparison we assume that the surveys cover 
most of the sky with \b\ > 45°. Weighting the Heiles model by sin \ b\ gives a prediction for what 
Lockman &: Gehman (1991) found; 

(JVgi”(|6|)sm|i.|) = JV“»sin(|(,|)df!. (42) 

= iVif‘f-l(l+sm|6|)JV£'r. (43) 

For |6| > 80°, LG91 find (A^g*)™(|6|) sin |6|) = 3.2 x lO^*’ cm“^ (note that the values they cite are 
for half the disk); the Heiles model gives 3.75 x 10^*^ cm“^, 17% larger. For |6| > 45°, LG91 find 
(Aig*)™(|6|) sin |6|) = 4.2 x 10^® cm~^, in excellent agreement with the prediction of the Heiles model, 
4.3 X 10^® cm“^. These results suggest that the Heiles model for optically thin HI is accurate to 
within about 20%. 

In a separate study, Lockman (1984) studied the HI in the inner Galaxy and concluded that 
three components were needed to describe the vertical structure, two Gaussians and one exponential; 
the exponential was required to account for gas high above the plane (see below). Kulkarni &: Heiles 
(1987) identified the component with the smallest scale height with the Cold Neutral Medium 
(CNM) and the other two components with the Warm Neutral Medium (WNM; we label these 
WNMi and WNM 2 ). In their review, Dickey and Lockman (1990) updated the properties of the 
three components. The total surface density in the Dickey-Lockman model, Aig*^™ = 6.2 x 10^® 
cm“^, is less than that directly measured outside the local “hole” by Heiles (1976), Aig*^™ = 7.45 x 
10^® cm“^. It is no surprise that these numbers are different: The Dickey-Lockman result applies to 
Galactocentric radii of 4-8 kpc, whereas the Heiles result is for the solar vicinity. Nonetheless, the 
two numbers agree to within 20%. More recently, Kalberla & Dedes (2008) found a mass surface 
density at the solar circle of = 10 Mq pc“^. Kalberla (private communication) has stated 

that this includes He, and that, furthermore, a more refined estimate gives = 9 Mq pc“^. 

This corresponds to A^g'^™ = 8.0 x 10^*^ cm“^, which is within 10% of the Heiles model. Since the 
Heiles model is intermediate between the Dickey-Lockman results and Kalberla’s, we shall adopt 
it as our fiducial model. 
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In order to infer the vertical distribution of the optically thin HI as well as the total column, 
we adopt the three components found by Dickey & Lockman (1990), but multiply the densities 
by 7.45/6.2 to bring the total column density up to the Heiles value. The central densities of 
the CNM and WNM in this model are 0.475 cm“^ and 0.206 cm“^; the effective scale heights, 
/leff = iV^*)™/(2no) are 113 pc and 338 pc, respectively. These values for the central densities are 
quite close to those found by Kalberla & Dedes (2008), 0.50 cm“^ and 0.19 cm“^, but the scale 
heights are somewhat smaller than they found, 150 pc and 400 pc. 

The total surface density of local HI in this model, = 6.0 Mq pc“^, is somewhat less than 
the maximum of about 10 Mq observed in star-forming disk galaxies (Wong &: Blitz 2002). This is 
reasonable, since the maximum is reached in regions that have significant amounts of molecular gas, 
whereas the local ratio of H 2 to optically thin HI is only about 0.12. The surface density including 
the He and heavier elements is 8.3 Mq pc“^. This is essentially the same as the 8Mq pc“^ given 
by Kulkarni &: Heiles (1987) after allowing for the difference in the He correction (they used 1.36 
instead of our 1.40); this is to be expected, since both values are based on Heiles (1976). We 
estimate the uncertainty in the surface density of the WNM as ~ 10% based on the magnitude of 
the stray radiation correction found by Heiles et al. (1981) at a typical point in the sky (\b\ ~ 30°). 
The CNM could have a larger error, since less than half the CNM within 1 kpc is in the volume 
with |6| > 10° with the density distribution in Table [2l we adopt 20%, including optical depth 
effects (see below). Since the columns of CNM and WNM are comparable and the errors could well 
be correlated, the total uncertainty in the HI column is about 15% 


5.2.1. Correction for Optical Depth 

The amount of HI given above is a lower limit since it does not allow for the opacity of the 
CNM, which causes the actual value of the HI column, N'm, to exceed the value inferred based 
on the assumption that the HI is optically thin, Strasser & Taylor (2004) found that 

TZ = AiHi/iVn*)™ = 1.32 based on the Canadian Galactic Plane Survey, which covered Galactic 
longitudes from 74° to 147° and Galactic latitudes from —3.6° to 5.6°. Braun (2012) has studied 
the HI in M31, M33 and the LMC and concluded that 77 = 1.34 ± 0.05 averaged over the entire 
galaxy. This value is in excellent agreement with that of Strasser &: Taylor (2004), although it 
must be borne in mind that the latter applies to gas close to the plane of the Galaxy, whereas 
Braun’s result applies to the entire galaxy. Braun (2012) attributed the self-absorption to large 
atomic clouds of order 100 pc in size. 

We can estimate the value of 77 for the solar vicinity from the results of Heiles & Troland 
(2003). Their data show that there is little absorption for \b\ > 40° (in fact, there is very little 
CNM for b > 40°), so we focus on the latitude range |6| < 40°. From Fig. 5 in their paper, we find 
that the 43 sources in this range have a mean value (77) = 1.29. They note that the line profiles 
for \b\ < 10° are complicated, which contributes to the uncertainty in this value. Unfortunately, 
there is no way to estimate the magnitude of this uncertainty; what we can conclude is that there 
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is no evidence that the local value of IZ differs significantly from the global values determined by 
Strasser &: Taylor (2004) and by Braun (2012). All these values are consistent with TZ = 1.30, and 
we adopt that value here. Note that TZ includes HI that is associated with partially molecular gas. 


Since it is only the CNM that is optically thick, we can determine the value of TZ in the CNM, 
"^CNM = ^CNm/-NcNM> 

'^CNmA'cNM + -^WNM 


7^ = 


^thir 

-'''hi 


In terms of the column density fraction /wNM = this implies 

^ — /wNM 


^CNM — 


1 — /wNM 


(44) 


(45) 


The Dickey-Lockman (1990) model has /wnm = 0.56; for TZ = 1.30, this implies that the CNM 
column densities are larger than the observed values by a factor T^cnm = 1-68. 

Our results for the local distribution of HI are summarized in Table [2l We have used T^cnm = 
1.68 to correct the observed values of the density and surface density of the CNM for optical depth 
effects. Since the three estimates we have for TZ are in good agreement, the uncertainty in TZ should 
not add significantly to the 20% error we have already estimated for the CNM. 


Recently, Fukui et al. (2015) have analyzed the all-sky LAB 21 cm survey (Kalberla et al. 
2005) together with the Planck map of of the dust emission at 353 GHz (Planck Collaboration 
2014) to infer the total amount of HI at |6| > 15°. A strength of their analysis is that the dust 
emission traces all the gas, including opaque HI. They reach two controversial conclusions: First, 
they infer that the “dark gas,” which gamma-ray observations show is not associated with either 
CO emission or 21 cm emission (Grenier et al 2005), consists primarily of opaque HI; by contrast, 
other observations (Lee et al. 2015 and references therein) and theory (Wolfire et al. 2010) are 
consistent with this gas consisting primarily of H 2 and atomic C. We note that the data they cite 
on direct measurements of H 2 absorption show that there is indeed little H 2 in regions in which 
A^hi ^ 1-0 X 10^^ cm~^; however, for the two cases in which the column density is larger, the H 2 
fraction is significant. Second, and of direct relevance to the present study, they conclude that 
typically 7^ ~ 2 — 2.5, suggesting that a significant amount of HI has been missed. However, this 
contradicts studies of HI absorption: Whereas Fukui et al. (2015) find a median value of the HI 
optical depth of order unity, Heiles & Troland (2003) find a median value less than 0.2 (Fig. 5 in 
their paper shows that the median value of i?raw = T/TZ exceeds 0.9, corresponding to thi < 0.21). 
In their study of the Perseus molecular cloud, Stanimirovic et al. (2014) find a median value of 
the maximum optical depth of thi = 0.16. The reason for this discrepancy is not clear. Reach et 
al. (2015) have confirmed the presence of significant amounts of dust that does not appear to be 
associated with CO or with HI emission in a set of interstellar clouds with masses of order 100 Mq. 
In order to account for the large amounts of dark gas in these clouds by optically thick HI, optical 
depths > 3 are required, which are rare. On the other hand, they show that a model in which the 
dust is associated with H 2 agrees with the the observed dust temperatures. They also discuss a 
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third possibility, that the observations are due to unusual dust rather than the presence of dark 
gas, and conclude that is unlikely. 

Regardless of the interpretation of Fukui et al.’s (2015) results, we find that they are in very 
good agreement with the Heiles model. Extrapolation of their results for the relation between dust 
emission and 21 cm emission for the warmest dust to zero 21-cm emission gives a negligible dust 
optical depth; we conclude that the ionized gas (Section l5.3p makes a negligible contribution to the 
dust. They excluded regions with obvious CO emission from their analysis, so the gas they analyze 
should consist of HI and dark H 2 that is not along the line of sight to CO emitting clouds. Locally, 
all the H 2 is only a small fraction of the HI, and if the H 2 along the line of sight to observable CO 
is excluded, the fraction is even smaller. We therefore expect gas observed by Fukui et al (2015) to 
be primarily HI. Under the assumption that the gas is located at a distance of 150 pc, their results 
correspond to a total mass of gas (including He) at \h\ > 15° of 1.4 x 10® Mq (Fukui 2015, personal 
communication). This corresponds to an average column density of HI as observed along a line of 
sight of (A^hi,1os)(|&| > 15°) = 6.0 X 10^® cm“^ (note that whereas A^hi is the column density through 
the entire disk, A^hi, los is measured from the Sun and goes through only half the disk). The average 
column density of optically thin HI in the Heiles model is (A^hi, ios)**^™(|('| > 15°) = 4.9 x 10^® cm~^. 
Since the WNM has a larger scale height than the CNM, gas with |6| > 15° has a lower fraction 
of CNM than the ISM as whole; hence, using TZ = 1.3 to correct for optical depth gives an upper 
limit on the average column density of HI, (A^hi, ios)(|5| > 15°) < 6.4 x 10^® cm“^. The average 
HI column density at |6| > 15° inferred from the Planck data by Fukui et al. (2015) thus agrees 
with the range of values in the Heiles model to within about 10%. This good agreement implies 
that our estimate for the mass of neutral gas (mostly HI in the interpretation of Fukui et al. 2015, 
but possibly including a significant amount of H 2 in other interpretations, such as that of Reach 
et al. 2015) in this part of the sky (excluding the small fraction with detectable CO emission) is 
reasonably accurate. However, the agreement does not explain why the Heiles model is dominated 
by optically thin HI whereas the analysis of Fukui et al. (2015) implies that a significant fraction 
of the HI is opaque. 


5.2.2. HI in the Galactic Halo 

For purposes of discussion, we define the gaseous Galactic halo in the solar vicinity to start 
1 kpc from the midplane. The Dickey-Lockman model gives A^hi(> 1 kpc) = 1.3 x 10^® cm“^; 
our renormalized version of this model gives 1.6 x 10^® cm“^. Lockman, Hobbs Sz Shull (1986) 
provided some support for the existence of HI far from the plane by comparing Lyman—a absorption 
observations of stars far from the plane with 21 cm emission observations in the same direction. 
They found that the HI column more than 1 kpc from the plane in the solar vicinity is (5 ± 3) x 
10^® cm“^. It shonld be noted that this result is less than 2 sigma, and that observations of 5 of the 
6 stars they observed more than 1 kpc from the plane were consistent with no HI in the halo. Their 
result is thus consistent with the lower values in the Dickey-Lockman and revised Dickey-Lockman 
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models. 

An independent estimate of the amount of HI in the local halo comes from 21-cm observations 
of high-velocity gas near the Galactic poles. Kulkarni & Fich (1985) found that up to 20% of the 
HI near the poles could be in a component with a velocity dispersion of 35 km s“^. Lockman 
&: Gehman (1991) modeled observations of the polar regions with hydrostatic equilibrium models 
and found that they needed a high velocity-dispersion component with a = 28, 34 km s“^ and an 
implied effective scale height fieflf = A'Hi/(2no) = 630, 770 pc for the two models they considered. 
If we approximate the vertical distribution of gas with an exponential of this scale height, then 
their results imply A^hi(> 1 kpc) ~ 1.4 x 10^® cm“^, where we have included the small contribution 
from the other components in their model. This result is thus consistent with both the original 
Dickey-Lockman model and our revision of it. 

Kalberla et al. (1998) found a larger amount of gas, with a total column density of 2.8 x 
10^® cm“^, with a higher velocity dispersion, a = 60 km s“^. They estimated the scale height of 
this gas as 4.4 kpc, so the amount more than 1 kpc from the plane is 2.2 x 10^® cm“^. This is 
somewhat larger than our estimate of 1.6 x 10^® cm“^, but the additional gas is more than 6 kpc 
from the plane. 

Our renormalized Dickey-Lockman model gives a column density of 1.25 x 10^® cm“^ above 
1.1 kpc, corresponding to a mass surface density of 0.14 Mq pc“®. This is a small correction, so 
the surface density of HI within 1.1 kpc of the plane is very nearly equal to the total value (Table 

ED. 


5.3. Ionized Gas: HII 

The total HII surface density at the solar circle can be separated into photoionized gas at ~ 10^ 
K and collisionally ionized gas at ~ 10® K (HIM). The photoionized gas can in turn be divided into 
diffuse ionized gas (DIG or WIM) and HII regions and their envelopes. Pulsar dispersion measures 
(DMs) provide a direct measurement of the column density of all the ionized gas. Since the HIM 
makes only a small contribution to the DM for z<5 kpc (Gaensler et al. 2008), we shall not treat it 
here. The classical model for the distribution of electrons in the Galaxy is that by Taylor & Cordes 
(1993), which was refined by Cordes & Lazio (2003). The dominant features in the solar vicinity 
in these models are a thick disk with a scale height ~ 0.9 kpc and the Gum Nebula. Subsequently, 
Gaensler et al. (2008) used pulsars with known distances to show that the scale height of the 
ionized gas is considerably larger, about 1.8 kpc. Recently, Schnitzeler (2012) has reviewed these 
and other models for the electron distribution and concluded that the model that best fits the 
observations is a revised Taylor-Cordes model in which the thick disk has a scale height of 1.59 
kpc and a DM normal to the plane of 24.4 cm“^ pc. Schnitzeler (2012) subtracted out the effects 
of all other features in the Taylor-Cordes model other than the thick disk before determining the 
parameters of the disk. The dominant such feature in the solar vicinity is the Gum Nebula, and 
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based on the work of Purcell et al. (2015), we estimate that this contributes about 0.11 Mq pc“^ 
to the surface density within a kpc of the Sun. We estimate that the mass of ionized gas in other 
features in the solar vicinity, such as the Eridanus Bubble, have much smaller masses of ionized 
gas. We present these results in Table [21 where we have added the 0.11 Mq pc“^ from the Gum 
Nebula to A^hIj ^ and i. An estimate of the error comes from comparing Schnitzeler’s DM with 
Gaensler et al.’s: they differ by only 5%. We have allowed for a somewhat larger fractional error 
in Si.i since the scale height is also uncertain. 


5.4. Local Interstellar Gas: Summary 

Table [2] summarizes our results for the three different components of the gas. Our results 
can be compared with the results of Bahcall et al. (1992), which were adopted with only minor 
changes by FOG. Bahcall et al. (1992) estimated 3 Mq pc“^ for the molecular gas based on the 
results of Scoville & Sanders (1987). These authors adopted a CO-to-H 2 conversion factor of 
3.6 X 10^^ cm“^ (K km s“^)“^; reducing that to the value recommended by Bolatto et al. (2013) 
reduces the column density of molecular gas to 1.7Mq pc“^. This is actually the value at the 
solar circle, not the solar vicinity, and since the Sun is in an interarm region it is not surprising 
that the local value we estimated, 1.0 Mq pc“^, is smaller. The HI data used by Bahcall et al. 
(1992) is based on Heiles (1976), just as ours is; we obtained 10.9 Mq pc“^ instead of 8 Mq pc“^ 
as they did because we corrected for optical depth effects and because we used a slightly larger He 
correction. We showed that a model for the distribution of local HI based on the work of Heiles and 
collaborators (the “Heiles model”) produces a total amount of local gas at |6| > 15° that is in good 
agreement with the results of Fukui et al. (2015) based on Planck data, implying that we are not 
missing a signihcant amount of unseen dark gas. For the ionized gas, Bahcall et al. (1992) adopted 
2 Mq pc“^ from Kulkarni & Heiles (1987), very close to the 1.8 Mq we obtained from more modern 
data. Since our somewhat larger value for the HI column compensates for our lower value of the 
H 2 column, we find a total surface density that is in good agreement with theirs: 13.7 Mq pc“^ 
vs. 13 Mq pc“^. However, this gas is often assumed to be entirely within 1.1 kpc of the plane, 
whereas more modern data on pulsar dispersion measures shows that almost half the ionized gas 
is more distant than that; as a result, we estimate that the surface density of gas within 1.1 kpc of 
the plane is 12.6 Mq pc“^. Finally, we note that Table [2] implies that the total midplane density 
of the local interstellar gas is reno = 1-17 cm~^, corresponding to pgo = 0.041 Mq pc“^, where we 
have summed the densities of each HI component rather than using the rounded value of 1.01 cm“^ 
to obtain this result. This is quite close to the value hno = 1-2 cm“^ obtained by Spitzer (1978), 
which was obtained by direct observations of gas in the plane. The uncertainty in this estimate is 
about 10%. The effective scale height of the gas is hg = Tig/{2pgQ) ~ 170 pc. 


Table 2. 


The Vertical Distribution of the Local Interstellar Medium'^ 



nniz) 

Nn,± 

y b 



(cm“^) 

(10^° cm“^) 

{Mq pc"2) 

(M© pc"2) 

H 2 

0.15 exp—(z/lOS pc)^ 

0.9 ±0.3 

1.0 ±0.3 

1.0 ±0.3 

HI: CNM'“ 

O.SOexp —(z/127 pc)^ 

5.54 

6.21 

6.21 

WNMi 

0.13exp —(2;/318 pc)^ 

2.24 

2.51 

2.51 

WNM 2 

0.077 exp —(z/403 pc) 

1.91 

2.14 

2.00 

Total HI 

?^HI,0 = 1-01 

9.7 ± 1.5 

10.9 ±1.6 

10.7 ± 1.6 

HH'i 

0.0154exp —( 2/1590 pc) 

1.6 ±0.1 

1.8 ±0.1 

0.9 ±0.1 

Total 

— 

12.2 ±1.5 

13.7 ± 1.6 

12.6 ±1.6 


“^References: H 2 : Dame et al. (1987), Heyer &: Dame (2015). HI: based on (see text) 
Heiles (1976), Dickey & Lockman (1990), Heiles &: Troland (2003). HII: Schnitzeler 
( 2012 ). 

^Includes He and heavier elements with 40% of the mass of H. 

'“Optically thin results for the CNM are q = 0.48 cm“^, ~ 3.30x 10^° cm“^, 

■^thin _ ■^thin _ o n/r 
^CNM ~ ^CNM,1.1 ~ ■ 

“^The values of Vh,_l, S, and Si.i include 0.11 Mq pc“^ for the Gum Nebula that is 
not included in hH(^:). 










- 28 - 


Table 3. The Local Baryon Budget {Mq pc ^), Excluding Halo Stars at z > 3 kpc 



S* 

E*i.i 


Egl.l 



Star Counts 

Elynn et al. (2006)*^ 

35.9^ 

32.0'^ 

13.2 

12.5'^ 

49.1 

44.5'^ 

Read (2014) 

37.2 ±1.2 

34.7 ± l.R 

17.0 ±4.4 

16.3 ±4.2^ 

54.2 ±4.9 

51.0 ±4.6® 

This work 

33.4 ± 3 

31.2 ±3 

13.7 ± 1.6 

12.6 ± 1.6 

47.1 ±3.4 

43.8 ±3.4 

Stellar Dynamics 

Kuijken k, Gilmore (1989b) 

35 ±5 

— 

13 ±3 

■ -- 

48 ±8 

— 

Garbari et al. (2012) 

— 

33.4iy 

— 

1 9+1.9 
-*-^-2.0 

— 

45.5+6 

Bovy k Rix (2013) 

38±4f 

— 


— 

51 ±4 

— 

Zhang et al. (2013) 

— 

42 ±5'^ 

— 

13® 

— 

55 ±5 

Bienayme et al. (2014) 

— 

— 

— 

— 

44.4 ±4.1 

43.3 ±4.0'^ 


'^Errors are taken from the original reference where possible, and do not always correspond to summing 
the errors in S* and in quadrature. 

’^Mixed approach; Stellar dynamics used to infer S* for main sequence stars with My < 8. 

■^Includes 0.4 Mq pc“^ for halo stars within 3 kpc of the plane. 

^Our inferred value based on data in reference. 

^Estimated value at 1.1 kpc based on our result, S*i.i = 0.934S*. 

^The stellar surface density is the total value, not the value within 1.1 kpc as implied by the text 
(Bovy, private communication). 

®Taken from E06; Zhang et al. (2013) assumed that this applies to gas within 1.1 kpc of the plane. 
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6. The Local Baryon Budget 

Our results for the stellar surface density, S*, gas surface density, and total baryonic mass 
surface density, Sb, in the solar vicinity are compared to those of others in Table [3l Total values 
(excluding stars at z > 3 kpc) and values within 1.1 kpc of the Galactic plane are given. We 
have included the results of Kuijken &: Gilmore (1989b), FOG, and all the “latest measurements” 
summarized by Read (2014), except for the following: Moni Bidin et al (2012b), whose work was 
criticized by Bovy &: Tremaine (2012); Bovy & Tremaine (2012), who investigated only the dark 
matter; and Smith et al. (2012), who deliberately presented an over-simplified toy model in which 
all the matter was divided into a thin sheet and a uniform distribution. We include the study 
of Bienayme et al. (2014), which appeared after the Read (2014) review. We have divided the 
references into two groups: The hrst group relies on star counts to infer the surface densities of 
the M dwarfs and stellar remnants (FOG, this work) or for all the stars (Read 2014). In the latter 
work, Read (2014) adopted = 30± 1 Mq pc“^ from Bovy et al. (2012b) and took the surface 

density of stellar remnants and brown dwarfs in the thin disk from FOG; he did not include the stellar 
remnants and brown dwarfs in the thick disk. The second group inferred the stellar surface density 
by using stellar dynamics to interpret observations of the stellar velocities of a tracer population. 

Our results for S* and Sb, both total and within 1.1 kpc, are consistent within 1 sigma with 
those of other workers with the exception of Read (2014) and Zhang et al. (2013). In inferring Sbi.i 
from Bienayme et al. (2014), we adopted their value for the scale height of 300 pc; had they used 
our results for the vertical distribution of baryons, the difference between their values of Sb and 
Sbi.i would have been larger. The discrepancy with Read (2014) is primarily due to the large gas 
surface density he adopted based on his assumption that the Kalberla &: Dedes (2008) HI column 
density did not include He, whereas it did (see above). The outlier among these works is Zhang 
et al. (2013), who find Sb = 55 ± 5 Mq pc“^. They provide only S*i,o, and we have made the 
conservative assumption that this is the same as S*i.i; if the vertical distribution of baryons follows 
our model instead, then their value for H*i,i would be 56 Mq pc“^. 

Summing the densities of the gas (Table [2]) and the stars (Table [T]) , we find that the total 
midplane density of baryons implied by our results is 

Pbo = 0.084 ± 0.012 Mq pc-3, (4G) 

where we have assumed an error of 15%, somewhat larger than that for the surface densities. This 
value of the midplane baryon density is somewhat less than (but within the errors of) the value of 
0.091 ± 0.009 Mq pc“^ found by FOG, primarily because their model of the interstellar gas gives a 
higher midplane density (hno = 1-45 cm“^ vs our value of 1.17 cm“^). Our value of pbo is within 
the errors of that of Bienayme et al. (2014), who found pbo = 0.077±0.007 Mq pc“^, and of Garbari 
et al. (2012), who found pbo = Mq pc“^. 
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7. The Local Density of Dark Matter, /3 dm 


It is now possible to combine our determination of the local surface density of baryons, 
with the measurements of the total mass surface density, S(z), by others to infer the local density 
of dark matter. There is no definitive evidence for a thin disk of dark matter (Kuijken &: Gilmore 
1989b), and existing data are consistent with an approximately spherical dark matter halo (Read 
2014). Here we therefore assume that the density of dark matter near the plane is constant, 
reserving discussion of a possible thin disk of dark matter to section \T?2\ below. It follows that the 
density of dark matter near the plane is proportional to the difference between the surface density 
of all matter and that of the baryons. 


^'DM 


S(z) - £/,(z) 
2z 


(47) 


The determination of S generally proceeds in two steps: First, the Boltzmann equation or an 
equation derived from it, such as the vertical Jeans equation, is used to analyze the kinematics of 
a tracer population of stars and thereby determine the vertical gravitational acceleration due to all 
matter within a distance z of the plane, K^- Then Poisson’s equation is used to relate S to Kz- As 
shown in Appendix [Cl this relation is 

\Kz 


S = 


27rG 


+ AS, 


where 


AS = 


1 d 


f 


v'^dz, 


(48) 
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27tGR or 

and where Vc is the circular velocity at {R,z). If the rotation curve is flat throughout the region 
of interest, then AS = 0 and Kz is directly proportional to S; this approximation is often made. 
However, Vc{R, z) generally varies with R for 2 7 ^ 0 even if the rotation curve is flat at z = 0 (Bovy 
&: Tremaine 2012 ). As a result, |Ar^| 7 ^ 27rGS even for a flat rotation curve. 

In Appendix Owe show that AS can be approximated as the sum of two terms. The first, 
ASi, depends on the slope of the rotation curve: 


ASi = 


azv: 


cO 


where 


a = 


d In Vco 


(50) 
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(9 Ini? 

measures the slope of the rotation curve and Vco is the rotation velocity in the midplane. This term 
vanishes for a flat rotation curve. Bovy et al. (2012a) find a = 0.01^q\, which is consistent with a 
flat rotation curve. The second term, AS 2 , is independent of the shape of the rotation curve and 
is generally ignored (e.g., Kuijken &: Gilmore 1989a). If the density in the disk (including the dark 
matter) is approximated as a double exponential, p oc exp(—i?//iij) exp(—z/Zi^), then this term is 


ASo = 


2 / 1 ^ 
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1 — e 


-z/h^ 


(52) 
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If the vertical distribution of matter were actually exponential, then Sq would be the total surface 
density of matter; this clearly does not work for the dark matter. We are interested in applying 
this equation at a distance ~ 1 kpc from the plane, which includes most of the baryons and where 
the surface density is dominated by the baryons. We therefore approximate Sq as the total surface 
density within a distance z ~ (1 — 1.1) kpc of the plane, and in our numerical evaluations we 
set So = 70 Mq pc“^. As noted in Appendix [Cl this double-exponential approximation for the 
density and surface density is accurate to within about 5%. Observe that AS 2 is positive and is 
independent of the rotation curve. We then have 

which clearly shows the dependence of S on the slope of the rotation curve. In AppendixOwe show 
that this result agrees well with the numerical results of Bovy &: Rix (2013). Now, R = 8.3 kpc 
(Chatzopoulos et al. 2015) and hji = 2.5 kpc (Bovy &: Rix 2013). We adopt the standard value 
for the rotational velocity at the solar circle, Vco = 220 km s“^, which is consistent with the recent 
result of Bovy et al. (2012a), Vco = 218 ± 6 km s“^. We borrow results from below to infer that 
hz = So/(2/Oo) — {70 Mq pc“^)/(2 x 0.093 Mq pc“^)=0.376 kpc. As a result, we obtain 

^ _ \Kz\ f 5.20ao.i + 2.00 Mq pc~^ {z = 1.0 kpc), 

2ttG y 5.72ao.i + 2.54 Mq pc“^ (z = 1.1 kpc), 

where ao.i = a/0.1. It should be borne in mind that both terms in AS increase with height above 
the plane: ASi oc z and AS 2 oc (approximately). 


7.1. Local Values of S and pdm 

Read (2014) has summarized the determinations of /?dm through 2013. Here we use our results 
for the baryon surface density, S^, and for the relation between S and Kz (Eq. (I54p i to revise the 
values of S and pdm for the solar neighborhood obtained by a subset of these references. We include 
the classic work of Kuijken k. Gilmore (1991). Among the latest local measures of pdm cited by 
Read (2014), we include Bovy k Rix (2013) and Zhang et al. (2013). As noted in Section El we 
exclude Moni Bidin et al. (2012b), Bovy k Tremaine (2012), and Smith et al. (2012). In addition, 
we have omitted Garbari et al. (2012), since Read (2014) has shown that the total surface density 
of matter inferred by their method is very sensitive to Two new determinations of pdm have 


^Garbari et al. (2012) find Eti.i = 45.5 M© pc“^ and pdm = 0.0221q q13 Mq pc“®, which implies a total column 
density of matter of Ei.i = 94 ± 31 Mq. On the other hand, Read (2014) shows that if one increases the surface 
density of baryons to 55 Mq pc“^, then poM drops from 0.022 Mq pc“® to 0.008 Mq pc“®, which corresponds to 
Ei.i = 74 Mq pc“^. In other words, increasing E;, by 10 Mq pc“^ decreases Ei.i by 20 Mq pc“^, and it reduces the 
implied poM by almost a factor 3. We emphasize that this in no way implies that their method is incorrect, only that 
we cannot use Equation ([47l) to infer the valne of poM for different values of E6 since E and E^ are strongly coupled. 
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appeared since Read’s (2014) review, one by Bienayme et al (2014), which we include, and one by 
Piffl et al. (2014). We do not include the Piffl et al. work in our analysis since their determination 
of Pdm is far more accurate (15% for a spherical halo) than an estimate based on Equation (1471) . 
We note that in many cases the distribution of matter is divided into a disk component, assumed 
to be baryonic, and a uniform component, assumed to be dark matter, and the masses of these 
components are determined from observations of stellar kinematics. However, as we have seen in 
Section [5^ the scale height of the ionized gas is sufficiently large (1590 pc) that it could be confused 
with the dark matter. Our method does not suffer from this potential problem since we determine 
the surface density of the baryons directly. 

Table m presents the data for the four studies we have selected to analyze to determine pdm- 
For Bovy &: Rix (2013), the value of S is taken directly from the reference; the other three references 
actually provide values of K^j (27rG), and we have converted to values of S by using Equation (f54]) 
under the assumption that the rotation curve is flat {a = 0). Inserting these values of S into 
Equation (|47p with our values of T,b, we obtain the values of pdm given in the right-hand column. 

Our results can be summarized as 

= 0.013 ± 0.003 M 0 pc-3 = 0.49 ±0.13 GeV cm’^, (55) 

where the uncertainty is the average of the uncertainties in Table 01 This uncertainty is based on 
the assumption that there is not a separate component of dark matter in a thin disk (see Section 
17.21 below). These uncertainties have been calculated from the quoted uncertainties in S given in 
the references together with the uncertainties in 'Eb from our work. One can show that the error 
estimate in Equation (|55p remains valid even if allowance is made for the uncertainty in the slope 
of the rotation curve found by Bovy et al. (2012a), —0.09 < a < 0.02. This result is consistent 
with that of Bienayme et al. (2014), pdm = 0.0143 Mq pc“^, but is more than a standard deviation 
above the estimates of Bovy & Rix (2013), pdm = 0.008Mq pc“^, and of Zhang et al. (2013), 
Pdm = 0.0065 Mq pc“^. The discrepancy with Bovy & Rix (2013) can be reduced to less than 1 
sigma by allowing for the fact that the column density of baryons within 1.1 kpc is less than the 
total; if the reduction is the same as we found, Ebi,i/Eb = 0.931, then the Bovy &: Rix value of 
Sfei.i would drop to 47.5 Mq pc“^, and pdm would increase to pdm = 0.0093 ± 0.0026 Mq pc“^. 
Comparing our result with recent references that are not included in Table 01 we note that our 
result is in excellent agreement with that of PifH et al. (2014), /?dm = 0.0126 ± 15% Mq pc“^, and 
it is consistent with the results of Garbari et al. (2012), pdm = 0-022lQg}3 Mq pc“^, and of Smith 
et al. ( 2012 ), pdm = 0.015 Mq pc“^. 

In our view, the cleanest determination of the local density of dark matter is that of Bovy &: 
Tremaine (2012, hereafter BT12), who analyzed the data of Moni Bidin et al. (2012a,b) on stars 
far enough above the plane (1-4 kpc) that the variation of gravitational acceleration with height 
is almost entirely due to dark matter. They found pdm = 0.008 ± 0.003 Mq pc“^. However, they 
then pointed out that there were two effects that would increase the density in the plane above 
this: First, since the dark matter was measured at an average distance of 2.5 kpc above the plane. 
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the radial distance to the Galactic Center is greater than in the midplane; in the NFW potential 
(Navarro et al. 1997) they adopted, this means that the density is 7% higher in the midplane 
than where it was measured. Second, in the gravitational field of the disk, the dark matter density 
varies as exp(—<l>/cj^), where <h ~ 2T:GT,hZ is the gravitational potential and a ~ 130 km s~^ is the 
velocity dispersion of the dark halo@ Altogether, this makes the local density of dark matter larger 
than that 2.5 kpc above the plane by a factor of about 1.3, so that 

Pdm(BT12) = 0.010 ± 0.004 Mq pc'^ (56) 

This is consistent with our result and with all the results in Table 01 

Read (2014) lists five global measures of the dark matter density based on fitting the rotation 
curve of the Galaxy under the assumption of a spherical dark matter halo. Our results are consistent 
with all of these, although only marginally in the case of Weber &: de Boer (2010), who found 
Pdm = (0.005 — 0.010) Mq pc“^. The two references that adopted the weakest priors are fully 
consistent with our result: Salucci et al. (2010) found pdm = 0.011 ± 0.005 Mq pc“^, and locco et 
al. (2011) found /?dm = (0.005 — 0.015) Mq pc“^. 

Finally, our estimate of the local matter density, both baryons and dark matter, in the plane 
is 

po = 0.097 ± 0.013 Mq pc-3. (57) 

This result is compared with others in Table [5l We include two classical pve-Hipparcos values, one 
by Oort (1960) and one by Kuijken &: Gilmore (1989); both are higher than our result. The highly 
cited work of Holmberg & Flynn (2000) has a value in good agreement with ours. Among Read’s 
(2014) “latest measurements,” only Garbari et al. (2012) give a result for pq; it is slightly more 
than one sigma above ours. Our result is in good agreement with the value obtained by Bienayme 
et al. (2014), pQ = 0.091 Mq pc“^. 


7.2. A Thin Disk of Dark Matter? 

Under the assumption that the density of dark matter near the plane is constant, we have found 
that the local ratio of the total density of matter to that of the baryons is po/pbo = 0.093/0.08 = 
1.16. Historically, some workers have found significantly larger values: for example, Oort (1960) 
found po/pbo = 2 (although stars fainter than My = 15 were not included in pbo); Bahcall et al. 
(1992) considered several models in which the dark matter was distributed like different components 
of matter in the disk and found po/pbo ^ 2 for the best fitting ones. 


®If the disk formed suddenly, then the velocity dispersion of the dark matter near the disk would be larger than 
this and there would be no overdensity. However, star formation appears to have proceeded steadily over the life of 
the disk (e.g., Binney et al. (2000), Just & Jahreifi 2010), so this effect is unlikely to be important. 
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Table 4. Local Density of Dark Matter 


Reference 

s. 


Pdm 

Pdm (This work) 


(Mq pC"2) 

(M© pc"2) 

(Mq pc"3) 

(Mq pc-3) 

Kuijken &: Gilmore (1991) 

Si.i = 73.5 ±6'^ 

Sfoi.i = 48 ± 8 

0.010® 

0.0135 ±0.0031 

Bovy & Rix (2013) 

Si.i = 68±4 

Sf, = 51 ± 4<^ 

0.008 

0.0110 ±0.0024 

Zhang et al. (2013) 

Si.o = 69 ± 6*’ 

Sw.o = 55 ± 5 

0.0065 

0.0130 ±0.0034 

Bienayme et al. (2014)® 

Si.o = 70.5 ± 1^ 

Sbi.o = 43 ± 4 

0.0143 

0.0137 ±0.0017 


®'From Equation (|47l) using the total surface density from the reference and the baryonic surface 
densities found here, Sfti i = 43.8 ± 3.4 Mq pc“^ and Sfti o = 43.1 ± 3.3 Mq pc“^. 

^Inferred from value of iL^/(27rG) given in reference using Equation (I54p with a = 0. 

'^Inferred from their values for Si.i and 

^The authors neglected the difference between and (Bovy, private communication). 
®Eor their scale height of 300 pc, Sbi.o — Sfei.i ~ 43 Mq pc“^. 


Table 5. The Oort Limit 


Reference po (-^0 pc 


Oort (lOOO)*" 

0.15 ±10% 

Kuijken Sz Gilmore (1989c)®'’^ 

0.11 - 0.29 

Holmberg &: Flynn (2000) 

0.102 ±0.010 

Garbari et al. (2012) 

U.iZU_o 019 

Bienayme et al. (2014) 

0.091 ±0.006 

This work 

0.097 ±0.013 


^Pie-Hipparcos estimates. 

^The intent of this work was to show that the estimation of the Oort limit depends on the star 
sample analyzed; the lower value corresponds to the estimate using a subsample of E8 dwarfs, 
whereas the upper one corresponds to a subsample of F5 dwarfs. 
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A high density of dark matter in the disk has been invoked to explain several terrestrial 
phenomena with periodicities of order 30 Myr. To our knowledge, the first such proposal was 
that of Rampino &: Stothers (1984), who suggested that the then-estimated periodicity in mass 
extinctions of 30 Myr (Raup &: Sepkoski 1984) could be naturally explained as the half-period 
of solar oscillations about the Galactic plane if po — 0.2 Mq pc“^, which was within the range 
of observed values at that time; they further suggested that “missing matter” might contribute 
to this density. Shaviv et al. (2014) inferred po = 0-21 Mq pc“^ in order to explain the 32 
Myr period observed in calcitic fossil shells; they attribute the excess density over the observed 
Po ~ 0.1 Mq pc“^ to a disk dark-matter component that is distinct from the standard halo dark 
matter. Fan et al. (2013) suggested that a more complex dark sector could include a small fraction 
of dark matter that is dissipative, so that it could cool and form a disk, like baryons do; they 
termed this “Double Disk Dark Matter.” Randall &: Reece (2014) explored a particular version 
of this model in which the disk of dissipative dark matter is very thin, so that comet showers 
are induced when the Sun passes through the disk with a half-period that matches the 35 Myr 
periodicity they infer from the record of large impact craters on Earth. In fact, as we shall show, 
all these models require that the disk of dark matter be significantly thinner than the baryonic 
disk. We define a “thin” dark disk as one with a thickness ^ that of the stellar disk (~ 400 pc, 
Eq. (I39l) l. This dark disk is therefore distinct from the one conjectured to be produced by accreted 
satellite galaxies, which could create a disk of dark matter with a density in the range 0.25 — 1.5 
times that of the halo dark matter but with a thickness > 1 kpc (Read et al. 2008, 2009). 

We obtain an approximate upper bound on the surface density of a possible disk of dark matter 
by comparing the dynamically determined disk surface density with that determined from star 
counts and the observed amount of gas (Kuijken & Gilmore 1989b, 1991). The determinations of 
Tif, in the “Stellar Dynamics” section of Table [3] are based on the assumption that the matter within 
a distance ~ 1 kpc of the Galactic plane is divided into a disk component, labeled “baryons,” and 
a dark matter component, which has a constant density. The “baryonic” component thus includes 
a possible disk of dark matter that has a scale height similar to or less than that of the baryons; 
we label the surface density of this component within a distance z of the midplane Sdisk, z- Our 
value for S;, is a valid measure of the baryons since it is based on (1) counts of the M dwarfs, (2) 
a determination of the surface density of the remaining stars that is consistent with the PDMF of 
RGH02, and (3) observations of the gas. Hence we can estimate the surface density of a possible 
disk of dark matter, SddMj by subtracting our determination of T,b from the dynamical ones in 
Table O Eddm = Ddisk, z — ^b,z- (Kramer &: Randall (2015) have shown that the presence of a 
thin dark disk can affect the determination of the stellar surface density; this effect is small for 
small values of SddM) and we ignore it for this approximate discussion.) Before carrying out this 
procedure, we note that the scale height of the ionized gas, h = 1590 pc, is so large that most of it 
would be counted as dark matter in a dynamical determination. Subtracting the column density of 
HII from the values we measure, we obtain S^^disk = 44.4Mq pc“^ and E;, disk,i.i = 42.1 Mq pc“^. 

The results of this exercise depend on the source of the data. The results of Bienayme et al. 
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(2014) imply that the surface density of disk dark matter is Sddm = 0 ± 5.3 Mq pc“^ (we have 
assumed that the disk surface density of 44.4 ± 4.1 Mq pc“^ that they found is the total for the 
disk, not the value for z = 1.1 kpc). At the other extreme, the high value of the disk surface density 
found by Zhang et al. (2013), 55.2 Mq pc“^ within 1.1 kpc, leads to Sddm = 13.1 ± 6.9 Mq pc“^. 
The average of all four estimates from the references in Table[3]is = 6.4±6.8Mq pc“^, which 

is consistent with both no thin disk of dark matter and with an average 1 sigma upper bound of 
Sddm < 13 Mq pc“^. This upper bound is adequate to accommodate the model of Randall and 
Reece (2014), which has Sddm — (10 — 13) Mq pc“^. 

In order to determine if a thin dark disk is consistent with stellar kinematic data, Kramer 
& Randall (2015) have repeated the analysis of Holmberg k. Flynn (2000) with updated data 
and an accounting for non-equilibrium effects. They point out that a thin disk of dark matter 
would gravitationally compress the stars, reducing the inferred surface density of stars with a given 
midplane density and velocity dispersion. Our determination of the stellar surface density is based 
on the midplane density and the observed scale height, so the sech^ part of the distribution is 
reduced somewhat less than they infer and the exponential part is unaffected. Furthermore, our 
determination of the surface density of gas is independent of the possible existence of such a dark 
disk. An accurate determination of the effect of a thin dark disk on the inferred surface density of 
baryons is beyond the scope of this paper, however. Kramer & Randall (2015) make the simplifying 
assumption that the baryons can be described by a single sech^ distribution. For an assumed scale 
height for the dissipative dark matter of 20 pc (they adopt a definition for the scale height that 
differs by a factor 2 from ours and give the scale height as 10 pc), they find a best-fit surface density 
of 4Mq pc“^ for the dark disk, which is less than the value needed to account for a terrestrial 30 
Myr periodicity. They also find a surface density of visible stars of 33 ± 2 Mq pc“^, marginally 
consistent with the result of Bovy et al. (2012b) but higher than the value we found above, and 
a total surface density within 1.1 kpc of 69 ± 6 Mq pc“^, in good agreement with the results in 
Table H 


As noted above, the scale height of a disk of dark matter with a midplane density comparable 
to the baryon density is significantly smaller than the baryonic scale height: 


^DDM = 


Sddm 


SpDM _ _ 

2/30,DDM 2(/3o,tot — PbO “ ^'halo Dm) ' 


(58) 


where /3o,tot includes the density of the disk dark matter and Phaio DM; which is denoted as /3 dm 
in the rest of the paper, is the midplane density of the halo dark matter. In order to bring the 
half-period of solar vertical oscillations down to about Pi 12 = 30 Myr, the total midplane density 
must be about 0.2 Mq pc“^ (Rampino k Stothers 1984) if the density is constant with height and 
greater than this if the scale height is not large compared to the amplitude of the solar motion. [§ 


®If the Sun has a vertical velocity of 7.25 km (Schonrich et al. 2010) and is 14 pc above the midplane (Binney et 
al 1997), then the amplitude of the solar oscillation about the midplane is 72[Pi/2/(30 Myr)] pc under the assumption 
that the motion is sinusoidal. 
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For Sddm < 13 Mq pc“^ this gives /iddm < 63 pc if we use our result (Eq. [55]1 for the density of 
the halo dark matter; in fact, the presence of disk dark matter would reduce the inferred density of 
the halo dark matter and slightly reduce /iddm- The model considered by Randall & Reece (2014), 
which has /iddm ~ (5 — 10) pc, readily satisfies this condition. A scale height < 60 pc would require 
a midplane density considerably greater than 0.2 Mq pc“^ in order to maintain a half-period of 
about 30 Myr; the numerical values of the Shaviv et al. (2014) model would require revision to 
accommodate the density variation in the dark disk. 

Before closing we note that while very thin dark disks are consistent with existing data, they 
face significant theoretical problems. The first is that the dark matter must have unusual properties, 
as in the model of Fan et al. (2013), in order to be able to collapse into a disk that is significantly 
thinner than the baryonic disk. Second, very thin sheets are cold and therefore gravitationally 
unstable (Toomre 1964). Galactic disks adjust themselves so that they are at least marginally 
stable against gravitational instabilities (e.g., Quirk 1972), so it is difficult for a self-gravitating 
structure to be significantly thinner, and therefore colder, than the gas in the disk. Further work 
on the physics of the proposed dark disks is required. 


8. Conclusions 

Stars. About half the stellar mass in the disk of the Milky Way is in the form of M dwarfs. 
The column density of M dwarfs in the solar neighborhood was determined by star counts in images 
taken with the Hubble Space Telescope (HST) in a series of papers by Bahcall, Flynn, Gould and 
collaborators, the last of which was Zheng et al. (2001; ZOl). We have revised the estimates of the 
surface density in those papers by using the determination of the density of stars in the immediate 
solar neighborhood (within 25 pc) by Reid et al. (2002; RGH02), obtaining an M-star surface 
density of 17.3 ± 2.3 Mq pc“^. 

We have determined the surface densities of main sequence stars more massive than M dwarfs 
by using the local densities determined by RGH02, the M-star scale height based on the results of 
ZOl, and the relative scale heights obtained by Flynn et al. (2006; F06). The overall normalization 
of the FOG scale heights appears too large since it disagrees with the values obtained by ZOl, Juric 
et al. (2008), Bochanski et al. (2010), Bovy et al. (2012b,c) and Bovy & Rix (2013). We therefore 
reduced the FOG scale heights by the factor required to bring the effective M-star scale height down 
to the ZOl value of 400 pc. This leads to a total surface density of visible stars (main sequence 
stars and giants) of 27 ± 2.7 Mq pc“^. This is marginally consistent with the result of Bovy et al. 
(2012b), who found 30 ± 1 Mq pc“^. However, the two results become fully consistent if allowance 
is made for the fact that their error estimate does not include the uncertainty in determining the 
mass range of the G stars they sampled. 

We have calculated the density of white dwarfs under the assumptions that the IMF above 
1 Mq is a power law, 'ip{m)dlnm oc m~^dlnm and that the star formation rate in the disk can be 


- 38 - 


approximated as a linear function of time. The white-dwarf density is relatively insensitive to the 
parameter describing this star formation history, bo, which is the ratio of the current rate to the 
average rate. For bo = 0.5, we find n*^wD = 8.5 x 10“^ pc“^. Katz et al. (2014) suggested that the 
density of white dwarfs has been significantly underestimated since the fraction of faint white dwarf 
companions to main sequence stars is much less than the fraction of bright white dwarf companions. 
Based on their discussion, we estimate that the density of white dwarfs is 8.3 x 10“^ pc~^, very 
close to the value obtained from direct calculation and larger than many previous estimates—e.g., 
6.0 X 10“^ pc“^ (Reid 2005) and 4.9 x 10“^ pc“^ (Sion et al. 2009). Since on average the progenitors 
of white dwarfs are older than a typical low-mass star, their scale height is somewhat larger (430 pc 
vs. 400 pc). For a mean white-dwarf mass of 0.665 Mq (Holberg et al. 2008), the surface density 
of white dwarfs is S*,wD = 4.9 Mq pc“^. 

Combining these estimates for the surface densities of main sequence stars, giants and white 
dwarfs with estimates of the surface densities of brown dwarfs, neutron stars and black holes (Table 
[I]), we find a total stellar surface density S* = 33.4 ± 3 Mq pc“^. This is close to the result of F06, 
who found S* = 35.9 Mq pc“^. Our results for visible stars are almost identical to theirs, although 
our mass distribution has significantly more M dwarfs than theirs; their estimate of the surface 
density of stellar remnants (mainly white dwarfs) is considerably larger than ours and therefore 
much larger than those of other workers, and they also estimated a surface density of brown dwarfs 
that is almost twice our value. Tabled] also gives the effective scale heights {h = T,/2po) for the 
different stellar and gas components, which can be useful in modeling the gravitational field normal 
to the plane in the solar vicinity. 

Gas. The properties of the local interstellar medium are summarized in Table [2] The principal 
uncertainty in the column density of interstellar gas has been the HI (Read 2014). We have 
carefully analyzed the existing observational data on the local HI and concluded that the local 
surface density of optically thin HI is 8.3 Mq pc“^ based on the observations of Heiles (1976) and 
Heiles et al. (1981), including He and heavier elements. This is essentially the same as the value 
of 8 Mq pc“^ given by Kulkarni Sz Heiles (1987) and is quite close to the value of 9 Mq pc“^ cited 
by Kalberla (private communication) based on the data in Kalberla & Dedes (2008). Based on the 
results of Heiles &: Troland (2003), Strasser &: Taylor (2004), and Braun (2012), we estimated that 
the optical depth in the 21 cm line implies that the total amount of HI is 1.3 times the optically 
thin value; for the CNM, which is the optically thick component of the HI, we inferred a correction 
factor of 1.68. We then found that the total surface density of local HI is 10.9 ± 1.6 Mq pc“^. This 
is significantly less than the value cited by Read (2014) of 17 ±4 Mq pc“^ since he took the surface 
density of HI to be 12 Mq pc“^ from Kalberla &: Dedes (2008) and, based on that reference, then 
added the mass of the He and heavier elements. We constructed a model for the distribution of 
local HI based on the work of Heiles and collaborators and showed that this model is consistent with 
the results inferred from Planck data by Fukui et al. (2015). This agreement has the important 
implication that our accounting of the local interstellar gas is not missing a significant amount of 
dark gas. Including the local surface density of H 2 and ionized hydrogen, we find a total gas surface 
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density of = 13.7 ± 1.6 Mq pc“^, in good agreement with the 13 Mq pc“^ given by Bahcall et 
al. (1992). The uncertainty in this value is significantly less than that in Read’s (2014) estimate. 
It is often assumed that all this gas is close to the Galactic plane, but the ionized hydrogen extends 
far from the plane (Gaensler et al. 2008, Schnitzler et al. 2012); we find that the surface density 
of gas within 1.1 kpc of the plane is 12.6 Mq pc“^. 

Baryons. Our results for the local surface densities of stars, gas and baryons are given in Table 
[3l together with the results of a number of other workers. Our result for the total baryon surface 
density, excluding halo stars at z > 3 kpc and hot halo gas, S;, = 47.1 ± 3.4 Mq pc“^, agrees 
with the results of the other references cited in the table with the exception of Read (2014), whose 
estimate of the gas surface density is much larger than ours as explained above. Our result for the 
baryon surface density within 1.1 kpc of the Galactic plane is significantly less than that of Zhang 
et al. (2012), who found S;, = 55 ± 5Mq pc“^. 

Dark Matter. The local mass density of dark matter can be determined by subtracting the 
surface density of baryons from the total surface density of matter determined from an analysis 
of stellar kinematics. In Table 0] we present the results from several references together with the 
results that would be obtained with our value for the baryon surface density. The average value 
we find from this exercise is /?dm = 0.013 ± 0.003Mq pc“^, which is consistent with the results 
of Bienayme et al. (2014) and Piffl et al. (2014), and with that of Bovy & Tremaine (2012) 
after including their corrections for the difference in the values of pdm in the midplane and at a 
distance of 2.5 kpc from the plane. Our result is also consistent with the global measures of the 
dark matter density discussed by Read (2014). We find that the total local density of matter is 
po = 0.097±0.013 Mq pc“^. The different determinations of the surface density of disk dark matter 
that we have analyzed lead to inconclusive results on the existence of a thin disk of dark matter. 
Using the method of Kuijken & Gilmore (1989b, 1991), we find that the average 1 sigma upper 
bound on a possible thin disk of dark matter in the Galactic plane is Bddm ^ 13 Mq pc“^, which 
is adequate to accommodate the models of Shaviv et al. (2014) and Randall & Reece (2014). This 
implies that if the midplane density is at least 0.2 Mq pc“^, which is necessary in order to match 
the the periodicities of several terrestrial phenomena, then the effective scale height of the disk dark 
matter is ^ 60 pc. This is consistent with the model of Randall &: Reece (2014) and of Kramer &: 
Randall (2015), but would require revision of the numerical parameters of the model of Shaviv et 
al. (2014). 

Determining the Surface Density S from the Aceeleration K^. Dynamical determinations of 
the total surface density of matter generally begin by determining the vertical acceleration due to 
gravity at some height z above the plane, K^. Appendix C gives the relation between S and 
and shows that it is in good agreement with the numerical results of Bovy & Rix (2013). 
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A. Effects of the Metallicity Gradient on Stellar Scale Heights 


The mean stellar metallicity decreases with distance from the midplane, which affects both the 
color-magnitude relation (CMR) used to estimate distances and the mass-luminosity relation used 
to infer stellar masses. Stars of a given luminosity are less massive at lower metallicity. Recall that 
we define the transition mass between M and K stars, tum-K) to be at My = 8, as in Gould et al 
(1996). This transition mass decreases with decreasing metallicity, dropping by about 17% at the 
point at which the metallicity is 0.5 dex below solar (Baraffe et al. 1998). 

Let z) be the number density of stars in the mass range m to m -|- dm at a distance 

z from the galactic plane; we denote the value at the plane as dn^,fi{m). For stars above the 
hydrogen-burning limit (m > mso)) we take m to be the Zero Age Main Sequence (ZAMS) mass, 
measured in solar masses. The number density of M dwarfs at height z is 




pmM- 

JiriBU 


mM-K{z) dn4rn,z) 


dlnm 


dlnm. 


(Al) 


The value at the midplane is and the corresponding mass density of M dwarfs at the midplane 

is p*,M,o- Let 


dN^,{m) dn^,{m,z) ^ 

dlnm dlnm 

'^max 


(A2) 


be the number surface density of stars per unit logarithmic mass interval within a distance |zmax| 
of the midplane. We shall assume that the stars are distributed symmetrically about the midplane. 
The number surface density of M dwarfs is then 


N, 


*,M 
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Jo Jrr 


2 / n^;udz, 
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^max dn^{m,z) 


dlnm 
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and the corresponding mass surface density is For stars with lifetimes greater than that of 

the disk, such as M dwarfs, we can express the number distribution in terms of the IMF (eq. [7]) as 


dn*(m) 


Fm j 


dlnm 


V^(m), 


(A5) 
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where 


Fm = 


rmM-K 
JrriBD 




(A6) 


is the number fraction of M dwarfs in the IMF at the midplane and where rriM-K = ?^m-k(0). 

The stellar metallicity gradient with distance from the plane implies that the vertical structure 
of the disk can be characterized by three distinct effective scale heights: 


• The effective scale height for a fixed mass, 


h{m) 


2dn^^o{m)/dliiim 


(A7) 


• The effective scale height for a fixed absolute magnitude; since we follow Zheng et al. (2001) 
in identifying stars with My > 8 as M dwarfs and since all M dwarfs have the same scale 
height, the M-dwarf scale height is 


hu = 


2n*,M,o ’ 


(A8) 


where is the number of M dwarfs per unit area identified by their absolute magnitudes. 


• The mass-weighted effective scale height, 

, _ 

where is the mass surface density of M dwarfs. 


(A9) 


One could also dehne an effective scale height based on the stellar colors, as done by Juric et al. 
(2008), for example. For a color-magnitude relation that is independent of height, as they assumed, 
this scale height corresponds to a fixed absolute magnitude and is similar to /im- In the following 
we show that the scale heights defined above are within a few percent of one another, so that the 
distinction among them, while important in principle, is not significant in practice. 

Given our assumption of a universal IMF, the effect of a variation in metallicity on the scale 
heights of M dwarfs is due to its effects on their mass range. The minimum mass of an M dwarf, 
ruBD) is very weakly dependent on metallicity (Burrows et al. 2001), so we neglect this variation. 
However, as noted above, the maximum mass of an M dwarf, itim-k, decreases by about 17% as 
the metallicity decreases from solar by 0.5 dex (Baraffe et al. 1998). For simplicity, we represent 
the density distribution of low-mass stars with a single exponential, 


With the aid of Equation 

A^*,m = 2n*,M,o^o — 


dn^{m,z) = dn*(m, 0) exp(—z/Zio) 

, Equation (IA4I) then becomes 
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Since the variation of mM-K(z) is not large, we obtain 


ho 


~ 1 — 


Fm/iq 
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^-z/ho 


mu-K 
_mM-K{z) 


dz. 


(A12) 


Using the results of Baraffe et al. (1998), ZOl modeled the metallicity gradient as a linear 
variation in the mass-luminosity relation between solar metallicity at the midplane and a metallicity 
lower by 0.5 dex at z = 1500 pc; they assumed a constant metallicity above 1500 pc. We implement 
their model by taking 


In 


rriM-K 

_mM-K{z) 
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rriM-K 


1500 pc [mM-K(1500 pc) 


(A13) 


We estimate that the logarithmic factor on the RHS is 0.16 from the results of Baraffe et al. (1998). 
As discussed in the text, ZOl found an effective scale height for M dwarfs of /im — 400 pc, and 
we adopt that value for ho- Since this is small compared with 1500 pc, we apply Equation (jA13|) 
beyond 1500 pc and let Zmax —>■ oo- Since tum-k = 0.67, Equation ()A12() becomes 


^ ^ _ 0.16V>(0.67) / ho \ 

ho Fm V 1500 pc/ 


Eor the standard PMHll IMF, ?/)(0.67) = 0.179 and Fm = 0.555, so that for ho = 400 pc. 


^ = 0.986. (A15) 
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One can repeat the same analysis for 
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and obtain hm,M/ho = S*^m/(2/?*,m,o) = 0.966. The ratio of hm,M to hyi, the quantity measured 
by ZOl, is 0.980. We conclude that these three different scale heights are approximately the same 
in practice. 


B. The Stellar Halo 


In their extensive study of the stellar number density distribution in the Milky Way based on 
data from the Sloan Digital Sky Survey, Juric et al. (2008) found clear evidence for a stellar halo. 
In our notation, their representation of the vertical stellar number-density distribution at the solar 
circle is 
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QhRo 
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{ (1 - I3)exp{-z/hi) + I3exp{-z/h2) + (1 -/3)A 
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The last term represents the stellar halo; the radius of the solar circle is Rq, which they took to 
be 8 kpc, fh is proportional to the fraction of stars in the halo, and the number measures the 
deviation from sphericity of the halo. The corresponding stellar surface density is 


Ai* = 2n 


*,0 


(1 - f3)hi + /3/i2 + -(1 - j3)fhqhRo'B 


- 

2 ’ 2 


(B2) 


where B(x,y) is the beta function. 

Juric et al. (2008) presented two hts to their data. Redder stars (r — i > 1.0) are too faint to 
have a significant halo contribution; we refer to the fit obtained from these data as the “disk ht,” 
which has //i = 0 by assumption. The scale heights are hi = 245 pc and /i 2 = 743 pc, and the 
thick/thin disk ratio is (3/{I — (3) = 0.13. The ht to the distribution of the bluer stars for R < 20 
kpc includes a halo component; we refer to this as the “disk + halo ht.” This ht has fh = 0.0051, 
qh = 0.64 (i.e., a battened halo), ly = 1.4, hi = 251 pc, /12 = 647 pc and (3/{l — (3) = 0.12. 
Corrections for binarity and Malmquist bias were presented only for the disk ht, so these values of 
the parameters do not include these corrections; it is not clear how these corrections would affect 
the inferred halo parameters. We assume that these corrections have at most a modest effect on 
the ratio of halo stars to disk stars. For example, halo stars contribute 8% of the total surface 
density in the disk + halo ht, and we anticipate that this result would be the very nearly the same 
if the corrected values of the scale heights, etc., were used. 

Juric et al. (2008) suggested that the disk ht includes many halo stars, and we have conhrmed 
that: the surface density for the two hts agrees to within 0.4% for 1 kpc < ^ < 3 kpc. At greater 
heights, the surface density of the disk + halo ht grows steadily larger than that of the disk 
ht, eventually being 4.4% larger. We conclude that including all the halo stars increases the 
surface density by about 4% above that obtained from the disk ht. This “high-halo” {z > 3 kpc) 
contribution includes about half the total number of halo stars. 


We assume that the halo was formed ~ 10 Gyr ago, so that only stars with My ^ 5 survive. 
The surface density of such stars in the disk is 23.1 Mq pc“^ according to Table 1, so the surface 
density of main sequence stars in the high halo is 4% of this, or 0.92 Mq pc“^. Brown dwarfs 
contribute 0.05 Mq pc“^. The contribution of white dwarfs is somewhat larger than 4% since all 
stars with 8>m>l have become white dwarfs; the parameter ai in Equation (jl7p is 0.7 for a 
starburst 10 Gyr ago instead of 0.53 for the case we considered in Section [3l We estimate a white- 
dwarf surface density of 0.35 Mq pc“^ for the high halo. The surface density of neutron stars and 
black holes in the high halo is negligible, so we ignore that. Altogether, we estimate that the high 
halo contributes a surface density of 1.3 Mq pc“^. The total halo contribution is about 2.6 Mq 
pc“^, but half of that is included in the disk hts. 

Most studies of the stellar surface density focus on stars below 3 kpc and assume that the 
density prohle can be described by one (e.g., Zhang et al. 2013) or two (e.g., Bochanski et al. 2010) 
exponentials, or the sum of an exponential and a sech^ function (e.g., ZOl). Halo stars within 
the survey volume are then automatically included in hts to the disk. We have therefore excluded 
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halo stars above 3 kpc from our analysis. By contrast, F06 explicitly included a halo component 
extending to the midplane. With an effective scale height h = S*/2/9*o = 3 kpc, about 2/3 of these 
halo stars are within 3 kpc of the disk, and these have been included in Table [TJ 


Although not of direct relevance to this paper, it is of interest to know the total mass of the 
stellar halo implied by the results of Juric et al. (2008). Their general expression for the stellar 
density in the halo is given by Equation (IB1|) with “1” in the final factor replaced by (R/Rq)'^, 
where R is the cylindrical radius. The mass in the stellar halo can then be obtained by integrating 
this equation over an ellipsoidal volume. Juric et al. (2008) state that their data extend out to 
galactocentric distances of about 20 kpc. Interpreting this distance as the height, Z, of an ellipsoidal 
volume, we find the mass of the halo within that volume to be 


Mh 


5,lxl0»(^)"(^)" 
\8 kpc J \20 kpc J 


3-2i/ 


Me 


0) 


(B3) 

(B4) 


for the parameter values adopted by Juric et al. (2008). In the final expression, we have set 
the midplane density of disk stars to be 0.0367Mq pc“^, which includes all the stars and stellar 
remnants with lifetimes > 10^^ yr from Table [TJ 


C. Relating the Total Surface Density, S, and the Gravitational Acceleration, 


Under the assumption of axisymmetry, Poisson’s equation, V-K = —AnGp, in cylindrical 
coordinates (i?, z) is 

dK, 1 d 2 


- = -dvrGp, 


(Cl) 


dz RdR 

where = —RKr is the circular velocity at R. Integration from —2 to 2 ; (assumed positive) gives 


K,- 


1 d r 
R^Jo 


vzdz = —27rGS 


(C2) 


(cf. Kuijken & Gilmore 1989a, who made the approximation that Vc is independent of z). Defining 


this becomes 
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1 d 
2ttGR'^ 
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2ttG 
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(C3) 

(C4) 


It is often stated that if the rotation curve is flat, then AS = 0 so that \Kz\ = 27rGS. However, 
this is not correct, since if the rotation curve is flat in the midplane, it is not flat away from the 
plane. Following Bovy & Tremaine (2012), we note that 

dKn d d(f) d d4> dK^ 
dz ^ ^ ^ 


(C5) 
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where cj) is the gravitational potential. Replacing Kr with —v'^jR and using Equation ()C4p . we 
obtain 


= VcO + 


(S - AS)dz', 


(C6) 


where Vco = Vc{z = 0) is the circular velocity in the midplane. Inserting this into Equation 
gives 


AS = 


1 d 


rz q pz' 

/ dz' vIq + 2txGR^ / (S-AS)(i 2 
4o Jo 


27rGR dR 

Generally, only the first term is retained (e.g., Kuijken &: Gilmore 1989a). 
Gonsider each of these terms in turn. The first term, labeled ASi, is 


(C7) 


ASi = 


azv, 
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where 


a = 


^GR?' 

d In Vco 


(C8) 


(C9) 


dlnR 

In their analysis of the Milky Way’s circular-velocity curve between 4 and 14 kpc, Bovy et al. 
(2012a) find 

(CIO) 


a = 0.0i;[|;?^ 


For a flat rotation curve, a and therefore ASi vanish, but the remaining two terms in Equation 
(|C7h do not, so that \Kz\ ^ 2 'kGTj. Numerically, ASi is less than 10% of the local surface density, 
S ~ 70 Mq pc“^, for z ^ 1.1 kpc: 


ASi = 5.72 


a0.l2l.l^c0,220 


Mq pc ^ 


RL 


(Cll) 


where ao.i = a/0.1, zi.i = ^;/(l.l kpc), Uco ,220 = Vco/i‘2'20 km s“^), etc. As a varies over the range 
allowed by observation, —5.15M0 pc“^ < ASi < 1.14 Mq pc“^ for nominal values of the other 
parameters. This term is often expressed in terms of the Oort constants A and B (Kuijken &: 
Gilmore 1989a), 

(A2 - B^)z 


ASi = -- 
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Next, consider the second term in Equation (IC7p . 
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1 d 52 
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RdR 9R2 
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We adopt a double-exponential form for the density, including the dark matter, in the solar neigh¬ 
borhood: 

p = p^^e^-R/hn-z/h.) ^ ^C14) 
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with hz constant; note that po = Poo exp(—The surface density is then 


S = Sooe"^/'*« (l - = Sq (l - , (C15) 

where Sq = 2po/i^ is the total surface density at R. Since the dark matter is not exponential, this 
dehnition has to be modified in practice. We are interested in applying this equation in the solar 
neighborhood at a distance ~ 1 kpc from the plane, which includes most of the baryons and where 
the surface density is dominated by the baryons. We therefore approximate So as the total surface 
density within a distance z ~ 1 — 1.1 kpc of the plane, and in our numerical evaluations we set 
Sq = 70Mq pc“^. With po = 0.093 Mq pc“^ (Eq. ISTll . the effective scale height is = 376 pc. 
Equation (jC13p then becomes 


ASs 



1 _ ^ ^ f 1 _ 


(C16) 


Note that AS 2 > 0 for the cases of greatest interest, in which R > and z > h^- For nominal 
values of the parameters {z = 1.1 kpc, hji = 2.5 kpc, R = 8.3 kpc and Sq = 70Mq pc“^), we 
find AS 2 = 2.54Mq pc“^. (We have found that applying the double-exponential model only to 
the baryons and using an NEW model for the dark matter changes this result by only about 5%.) 
Hence, we have 

— 2.6 Mq pc“^ < (ASi -|- AS2) < 3.7 Mq pc“^ (C17) 


after allowing for the observed range of o. (Eq. IClOp . 


Finally, we show that the third term in Equation (IC7h . 


AS 3 


\RdR ^ 
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is negligible. This is to be expected, since a rough self-consistency argument shows that AS 3 is 
second order in AS/S: 



(C19) 


where the first step follows from Equation (IC7p and the second step follows from the fact that AS 2 
is of order AS. Because AS = ASi -|- AS 2 -|- AS 3 , Equation (jClSp becomes 


AS 3 — AS 31 -|- AS 32 -I- AS 33 


(C20) 


in an obvious notation. Since the rotation curve is approximately flat, we neglect the radial variation 
of the velocity and find 


AS31 ~ --(1 - af 



ASi ~ -O.OI2AS1, 


(C21) 


which is in the range —0.014 Mq pc AS 31 < 0.06 Mq pc ^ for nominal values of the parameters 
and is negligible. An upper limit on AS 2 is 


AS 2 < 


(z^ ^ 2hl\ 


So 


(C22) 
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from Equation (jC16l) . We then obtain 


AS32I < 


/ + 12hlz‘^\ 

V 244 J 


So. 


(C23) 


For nominal values of the parameters, this upper bound is 0.26 Mq pc which again is negligible. 


To evaluate AS 33 , we must specify the radial and vertical dependence of the unknown term 
AS 3 . We expect the radial dependence AS 3 to be qualitatively similar to that of ASi, which 
varies on the scale of the half-radius, R/2 (since ASi oc R~^), and AS 2 , which varies on the radial 
scale length, Hr. Since Hr ~ 2.5 kpc < R/2 ~ 4 kpc, we obtain an approximate upper bound on 
AS 3 by evaluating the radial derivatives with the radial scale length. Since AS vanishes for z = 0 
according to Equation (|C3I) . we assume that AS 3 oc z'^ for some value of p > 0. Equation (jClSp 
then implies 


AS 


22AS3 


33 


{p + 1)(7' + 

Since AS 33 has opposite sign to AS 3 , Equation ()C20p implies that 


(C24) 


IAS 3 I < IAS 31 -|-AS 32 I. (C25) 

For nominal values of the parameters, we find IAS 3 I <O.27M0 pc“^, which is negligible. 

We conclude that 

S~S + ASi + AS 2 , (C26) 

ZTTCt 

where —2.6 Mq pc“^ < (ASi -|- AS 2 ) < 3.7 Mq pc“^ (Eq. IC17P for fiducial values of the parame¬ 
ters as a varies over the allowed range (Eq. IClOp . 

We can test our result by comparing with the results of Bovy & Rix (2013). We begin with 
their fiducial case, in which the rotation curve is flat locally (a = 0). They give the values of S and 
Kz for a number of mono-abundance populations measured at different Galactocentric radii. For 
the bin with data that is closest to the Sun (7.5 kpc < R < 8 kpc, where they adopted Rq = 8 kpc), 
the average value of S — |it'^|/(27rG) is 2.6 Mq pc“^. Using their values of (400 pc) and Hr (2.5 
kpc, which we have adopted as our fiducial value) and using the average value of Eq in this bin 
(73.4 Mq pc“^), we find AS ~ AS 2 = 2.5Mq pc“^ from Equation ()C16I1 . in excellent agreement 
with their result. In making this comparison, we have ignored the small correction due to the fact 
that the slope of the rotation curve varies with radius in the Bovy & Rix (2013) model (Jo Bovy, 
private communication). 

They also considered the effects of non-flat rotation curves with a = ±0.1 at the solar circle. 
They adopted Uco = 230 km s“^, so that ASi = ±6.7Mq pc“^ for a = ±0.1. For a = —0.1, 
their analytic fits give Ur = 2.3 kpc and S = 65 Mq pc“^ and |R^|/(27rG) = 69 Mq pc“^, 
which implies AS 2 = 2.75 Mq pc“^. Given that S = 65 Mq pc“^, we obtain |R'2|/(27rG) = 
69.0 Mq pc“^, in excellent agreement with their result. For a = ±0.1 they found Hr = 2.7 kpc 
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and S = 72 Mq pc“^, |i^^|/(27rG) = 63 Mq pc“^. (The difference between |ii'2|/(27rG) and S is 
larger than given in Equation (|C17I) primarily since +0.1 is outside the observed range of a given 
by Bovy et al. (2012a).) We then obtain AS 2 = 2.06Mq pc“^; for S = 72 Mq pc“^, our results 
give \Kz\/{27:0) = 63.2Mq pc“^, which again is in excellent agreement with their result. 
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